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ABSTRACT 


An  investigation  is  being  conducted  to  determine  whether 
ballistic  damage  can  seriously  degrade  the  aeroelastic  integrity  of 
lifting  surfaces  on  aircraft.  A  potential  aeroelastic  failure  mode 
that  was  identified  in  the  first  year's  study  has  been  investigated 
here  over  a  larger  range  of  parameters.  This  failure  mechanism  results 
from  the  localized  steady  drag  generated  when  a  lifting  surface  encounters 
damage  to  its  aerodynamic  shape.  Its  modeling  has  been  extended  in 
this  study  to  swept  wing  configurations  and  to  possible  multiple  and 
distributed  damage  sites.  In  addition,  a  larger  range  of  single  damage 
site  locations  have  also  been  considered  to  assess  the  possible  trade¬ 
offs  between  the  influence  of  both  structural  and  aerodynamic  damage 
locations.  A  check  on  the  validity  of  the  strip  theory  aerodynamic 
modelling  employed  in  this  study  has  also  been  made  by  comparing  these 
results  with  those  obtained  from  a  lifting  surface  theory  modeling. 
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Finally,  an  additional  failure  mechanism  is  identified  that  results 
from  any  unsteady  but  periodic  fluctuating  aerodynamic  drag  loads  that 
are  generated  by  the  damage.  A  parametric  and  oscillatory  instability 
can  be  induced  by  relatively  low  level  drag  loads  in  this  case  if  they 
happened  to  be  appropriately  tuned  to  the  structural  frequencies  of  the 
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1 .  INTRODUCTION 


The  influence  of  ballistic  damage  on  the  aeroelastic  response 
of  an  aircraft  in  flight  is  not  yet  well  understood.  The  primary 
damages  of  interest  are  those  sustained  by  the  lifting  surfaces,  that 
is,  the  aircraft  wings  and  tail  sections.  A  review  of  inflight  films 
illustrating  ballistic  damage  to  aircraft  indicates  that  in  some  cases 
an  aircraft  can  tolerate  a  significant  amount  of  damage  from  several 
hits  while  in  other  cases  a  single  hit  may  result  in  the  immediate 
destruction  of  the  aircraft.  A  question  naturally  arises  as  to  whether 
a  reasonably  small  amount  of  damage  might  occur  in  a  critical  area  that 
could  promote  an  explosive  type  of  flutter  or  divergence  instability. 
These  instabilities  would  most  likely  destroy  the  aircraft.  In  the 
present  study  an  investigation  is  carried  out  to  determine  whether 
ballistic  damage  can  seriously  degrade  the  aeroelastic  integrity  of 
lifting  surfaces  on  aircraft. 

Some  of  the  first  and  most  comprehensive  work  on  this  general 
subject  was  conducted  in  1950  by  Biot  and  Arnold  as  out-lined  in 
Reference  1.  The  results  of  their  studies  demonstrated  that  aeroelastic 
instabilities  were  not  easily  triggered  by  ballistic  damage.  Further¬ 
more,  if  sufficient  structural  damage  was  imposed  on  a  lifting  surface 
to  lower  its  flutter  and  divergence  speeds  into  the  flight  envelope, 
the  surfaces  would  fail  due  to  inadequate  strength  rather  than  due  to 
in  adequate  stiffness.  In  essence,  the  reduction  of  flutter  and/or 
divergence  margins  by  25%  required  nearly  an  90%  loss  of  torsional  stiff- 
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ness  at  a  certain  critical  section  within  the  wing.  It  was  probably  the 
impact  of  this  finding  that  delayed  any  further  investigation  of  this 
subject  for  nearly  30  years. 

A  recent  investigation  by  Hemmig,  Venkayya,  and  Eastep  (Reference 
2)  has  incorporated  more  contemporary  finite  element  techniques  to  model 
structural  damage  to  lifting  surfaces.  Their  results  also  suggest  the 
difficulty  of  reducing  flutter  and  divergence  margins  of  highly  redundant 
lifting  surface  structures  much  below  the  levels  suggested  by  Biot  and 
Arnold.  Some  further  insight  into  the  problem  is  obtained  here,  however, 
when  it  is  demonstrated  that  highly  optimized  structural  designs  may 
demonstrate  increased  aeroelastic  sensitivity  to  ballistic  damage. 

Recent  aeroelastic  investigations  of  lifting  surfaces  within  the 
literature  suggests  that  refinements  in  the  aerodynamic  modeling  may  be 
in  order  (3,4).  In  essence,  the  chordwise  forces  due  to  drag  and  leading 
edge  suction  are  found  to  have  an  observable  influence  on  the  flutter 
and  divergence  boundaries  for  certain  wing  geometries.  In  addition, 
other  aerodynamic  investigations  outlined  in  References  5  and  6  suggest 
that  wings  with  a  highly  optimized  aerodynamic  configuration  may  demon¬ 
strate  increased  sensitivity  to  ballistic  damage  through  larger  drag 
rises.  In  view  of  this,  the  present  study  concentrates  on  an  extension 
of  the  earlier  Biot  and  Arnold  work  by  incorporating  into  their  analysis 
chordwise  forces  due  to  drag  that  arise  as  the  result  of  the  ballistic 
damage. 

A  potential  aeroelastic  failure  mechanism  was  identified  in  the 
first  year’s  study  that  results  from  the  localized  drag  that  can  be 
generated  on  a  lifting  surface  due  to  significant  damage  to  its  aero¬ 
dynamic  shape.  When  this  localized  drag  occurs  at  critical  positions 
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over  the  surface,  it  drastically  lowers  the  divergence  speed  of  the  sur¬ 
face.  For  this  reason,  this  failure  mechanism  will  be  referred  to  as 
a  "drag  divergence"  mode  of  instability.  The  present  study  extends 
this  drag  divergence  investigation  to  a  larger  range  of  parameters.  This 
includes  several  damage  site  locations  to  assess  the  possible  trade¬ 
offs  between  the  influence  of  both  structural  and  aerodynamic  damage 
locations  on  the  wing's  divergence  characteristics. 

A  second  study  is  made  to  determine  the  influence  of  ballistic 
damage  on  more  contemporary  swept  wing  configurations.  Both  swept  back 
and  swept  forward  wings  are  considered  that  have  the  same  stiffness 
distribution  along  the  elastic  axis  as  the  statistical  straight  wing  model 
studied  earlier. 

The  previous  drag  divergence  studies  were  also  extended  to  allow 
for  the  possibility  of  distributed  drag  effects  and  for  multiple  con¬ 
centration  drag  influences  on  both  straight  and  swept  wings.  This 
extension  allows,  for  example,  the  inclusion  of  the  distributed  un¬ 
damaged  wing  drag,  store  drag,  gun  cannon  recoil,  and  multiple  damage 
site  modeling  in  the  aeroelastic  analysis. 

A  second  aeroelastic  failure  mechanism  is  identified  for  the 
case  when  the  damage  gives  rise  to  periodic  aerodynamic  forces  in  the 
chordwise  plane  of  the  lifting  surface.  When  these  periodic  forces 
are  appropriately  tuned  to  the  structural  mode  frequencies  they  can 
force  the  wing  into  a  dynamically  unstable  oscillation  known  as  parametric 
resonance.  This  mode  of  instability  is  demonstrated  on  a  statistical 
model  of  a  fighter  wing. 

Finally,  it  should  be  mentioned  that  the  transient  structural 
responses  of  the  lifting  surfaces  to  the  pressure  forces  arising  from 
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the  explosive  impact  of  the  ballistic  warhead  are  not  considered  in  the 
present  drag  divergence  study.  In  essence,  only  those  structural  responses 
are  investigated  that  are  due  to  the  steady  state  or  periodic  air  loads 
imposed  by  the  damage  and  the  assumed  harmonic  unsteady  or  quasi-steady 
air  loads  that  are  classically  imposed  in  a  standard  flutter  analysis. 


2.  DAMAGE  IDENTIFICATION 


Two  types  of  aircraft  lifting  surface  damage  are  considered 
in  this  report.  One  is  the  structural  damage  which  is  reflected  as 
a  reduction  of  the  bending  stiffness  (El)  and  torsional  stiffness  (GJ) 
of  the  lifting  surfaces.  The  other  is  aerodynamic  damage  which  is 
the  modification  of  the  aerodynamic  forces  arising  from  the  change  of 
shape  of  the  aircraft  lifting  surfaces. 

2.1  Structural  Damage  Model! ing 

In  developing  a  model  for  structural  damage  studies,  reference 
at  this  point  is  made  to  the  comprehensive  work  by  Biot  and  Arnold 
(Reference  1),  which  assesses  the  changes  in  the  physical  parameters 
of  a  wing,  such  as  rigidity  changes,  and  their  secondary  effect  on 
inertia  due  to  a  shift  in  elastic  axis  that  might  be  caused  by  ballistic 
damage.  In  that  study,  the  authors  conduct  an  investigation  which  in¬ 
volves  the  flutter  analysis  of  a  typical  model  of  a  fighter  wing  subject 
to  a  wide  range  of  structural  damage  conditions.  The  basic  aeroelastic 
parameters  of  this  typical  wing  are  chosen  to  be  the  geometric  mean  of 
some  fourteen  different  fighter  aircraft.  These  are  representative  of 
early  1950  fighter  configurations  which  have  moderate  aspect  ratio 
straight  wings.  It  was  further  de'  onstrated  in  this  study  by  Biot  that 
the  statistical  model  reflected  the  correct  trends  in  modification  of 
aeroelastic  characteristics  due  to  damage  that  was  determined  for  a 
randomly  selected  fighter  from  the  group.  Table  1  contains  a  listing 
of  these  fourteen  aircraft.  The  geometric  mean  of  the  aeroelastic 


parameters  of  this  typical  model  are  given  in  Table  2.  Since  aircraft 
of  this  type  are  normally  designed  to  withstand  loads  up  to  1.5  times 
the  maximum  limit  load  (load  factor  3.0  to  8.67),  the  structure 
should  normally  experience  decreases  in  torsion  rigidity  (GJ) 
and  bending  stiffness  (El)  of  up  to  45%  and  not  undergo  strength 
failure  at  its  limit  load  factor.  Higher  percentage  could,  of 
course,  be  tolerated  at  lower  load  factor.  To  model  in  detail 
the  changes  ir.  structural  parameters  due  to  damage  the  two  spar 
wing  is  idealized  as  illustrated  below. 


In  the  analysis  of  such  a  two  spar  wing  it  is  common  practice  to 
neglect  all  .  f  the  material  aft  of  the  rear  spar  where  the  controls 
are  generally  located.  A  parametric  study  is  conducted  on  such  a 
model  in  Reference  1  to  determine  the  shift  that  could  be  imposed 
on  the  wing  elastic  axis  by  ballistic  damage.  Such  a  modeling 
indicates,  for  example,  that  even  extreme  damage  patterns  result 
in  a  movement  of  the  elastic  center  of  the  wing  by  only  12.5%  of 


the  section  semi-chord.  In  essence,  it  is  found  that  shifts  in  the  elas¬ 
tic  axis  due  to  ballistic  damage,  which  still  left  the  structure  with 
adequate  strength,  would  not  significantly  influence  the  flutter  and 
divergence  margins.  Consequently,  in  the  structural  damage  model  em¬ 
ployed  for  the  present  study  only  reductions  of  El  and  GJ  are  applied  in 
the  damage  area  to  reflect  ballistic  damage  to  the  structure.  Since  the 
wing  structure  of  this  typical  fighter  is  a  two  spar  semi-monocoque 
construction  (see  Figure  1),  torque  box  destruction  and  GJ  reduction 
would  probably  be  easier  to  accomplish  than  reducing  El  by  removing 
spar  cap  material.  For  this  reason,  the  more  representive  structural 
damage  cases  are  thought  to  be  those  involving  larger  GJ  reduction  than 
El  reductions. 

It  is  apparent  that  the  finite  element  method  can  be  employed 
in  a  more  detailed  structural  modelling  when  a  smaller  aspect  ratio  wing 
must  be  identified  or  when  more  structural  detail  is  available  for  a 
statistical  model.  The  reader  is  referred  to  Appendix  C  and  Reference  2 
for  more  information  on  this  type  of  modelling. 

2.2  Aerodynamic  Damage  Modelling 

The  present  investigation  employs  basically  the  same  strip 
theory  aerodynamic  modelling  as  utilized  in  Reference  1.  However,  one 
additional  parameter  (not  in  Reference  1)  is  included  in  this  study. 

This  parameter  accounts  for  the  steady  state  chordwise  concentrated  drag 
force  that  occurs  as  a  result  of  aerodynamic  damage.  The  influence  of 


possible  time  dependent  unsteady  dt ag  effects  are  discussed  in  Section 
7  of  this  report.  This  first  order  influence  of  the  aerodynamic  damage 
is  based  upon  an  assumed  localized  hole  through  the  lifting  surface. 


Furthermore,  this  damage  is  considered  to  introduce  only  local  or  con¬ 
centrated  changes  in  the  aerodynamic  lift,  drag  and  moment  distributions 
as  illustrated  in  Figure  2.  This  localized  modelling  of  the  loads 
induced  by  the  thru  hole  type  damage  is  justified  based  upon  the  find¬ 
ings  from  a  kernel  function  lifting  surface  theory  study  of  these  types 
of  discontinuities  (Reference  7).  Figure  3  illustrates  such  an  example 
taken  from  a  theoretical  study  in  Reference  7  where  a  thru  hole  type 
damage,  shaded  region,  produces  chordwise  and  spanwise  pressure  distur¬ 
bances  that  die  out  within  a  characteristic  hole  dimension.  The  experi¬ 
mental  observations  in  Volume  II  of  this  report  also  confirm  these  findi 
A  coupling  of  these  damaged  induced  loads  with  the  thin  beam  lateral 
buckling  equations  for  the  wing  structure  results  in  only  the  drag  pro¬ 
viding  the  first  order  aeroelastic  effect  of  the  damage.  It  was  anti¬ 
cipated  that  this  concentrated  type  of  aerodynamic  drag  force  might  be 
capable  of  producing  a  structural  failure  of  the  wing  even  in  the 
absence  of  significant  structural  damage.  This  is  illustrated  in  Figure 
4  which  indicates  that  high  chordwise  drag  loads  generated  near  the  tip 
of  the  wing  may  cause  the  wing  to  snap  or  diverge  laterally  similar  to 
the  lateral  buckling  of  a  thin  beam. 

An  estimate  of  potential  damage-induced  drag  levels  can  be  ob¬ 
tained  from  the  experimental  studies  of  Reference  8  or  9,  where  changes 
due  to  damage  in  the  lifting  surface  drag  polar 

CD  =  Cq  +  PCL 

TT 

are  expressed  as  changes  in  the  parasite  drag  coefficient 

TT 

and  the  slope  p  of  the  induced  drag  term.  A  coefficient  of  drag 
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increase  is  defined  in  terms  of  a  damage  area  as 

S  =  wing  area 
A  =  hole  or  damage  area 


That  is,  CD  and  p  represent  increments  in  the  drag  polar  based 

upon  the  hole  damage  area,  i.e.,  the  sum  of  the  entrance  and  exit 

hole  areas  divided  by  two.  A  tabulation  of  experimentally  deter- 
★  ★ 

mined  values  of  and  p  is  given  in  Table  3.  These  parameters 
are  measured  on  a  2-ft.  chord  and  an  8^-ft.  span  two-dimensional 
wing  section  having  a  symmetrical  NACA  65^ -012  profile.  The 
Reynolds  number  of  the  tests  is  3.7  x  10^  which  is  sufficiently 
large  that  the  measurements  may  be  applied  to  full-scale  damaged 
aircraft.  The  classes  of  damage  (i.e.,  leading  edge,  midchord, 
etc.)  are  also  presented  in  Table  3.  From  this  table  it  is 
evident  that 


*  < 

p  16 


represent  experimentally  determined  upper  bound  values  for  the 
drag  parameters.  These  larger  values  occur  for  the  leading  edge 


and  then 


with 


r*  =  AD  _  ,r  S 
"  qA  "  A 


*  *  .  ; 
CD  =  CD  +  P  C[ 


*  8 
C  =  C  — 

ld  a 

IT  IT 


P*  =  AP  ^ 


class  A  or  B  type  damage.  Based  upon  the  data  of  Table  3,  a 
conservative  damage-induced  drag  estimate  would  be 
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1.0  +  5C* 


while  an  upper  bound  estimate  for  the  higher  Mach  numbers  and 
larger  hole  sizes  would  appear  to  be 

C*  =  2.4  +  10C^ 


These  damaged-induced  drag  coefficients  are  based  upon  the 
assumption  that  leading  edge  type  damages  are  imposed  upon  the 
lifting  surfaces. 

To  obtain  an  estimate  of  the  actual  levels  of  damage- 
induced  drag  force  *hat  can  be  imposed  on  the  typical  fighter 
aircraft  in  combat,  consider  the  following  fighter  parameters. 


weight  w  =  21,000  lb  (including  stores) 

2 

wing  area  S  =  300  ft  (including  fuselage  carry  thru) 

aspect  ratio  =  4.5 

Cp  =  0.014  (clean  fighter) 

utt 

Cp  =  0.026  (fighter  with  external  stores) 

7T 

p  =  0.083 


combat  altitude  20,000  ft 


Based  upon  these  parameters,  an  assumed  600  mph  flight  speed, 
and  the  experiments  of  Reference  9,  one  can  estimate  the  range 

i 

I 

of  drag  increments  that  would  be  expected  for  a  range  of  load 
factor  'n‘  and  damaged  area  ‘A1  expressed  in  terms  of 
%  wing  area  for  several  combat  altitudes.  These  results  are 
presented  in  Table  4  for  20,000  ft  altitude  and  in  Figures  5 
and  6  for  the  different  combat  altitudes. 

Several  other  drag  factors  that  are  not  considered 

here  could  also  be  simultaneously  degrading  the  divergence 
characteristics  of  the  wing.  External  stores  on  the  wing,  for 
example,  give  rise  to  concentrated  drag  and  chordwise  loadings 
as  does  the  recoil  from  firing  wing  mounted  gun  systems.  The 
drag  on  fighters  loaded  with  external  stores  can  be  nearly 
double  the  drag  of  a  clean  fighter  in  a  cruise  condition.  In 
addition,  the  recoil  loads  on  a  wing  due  to  rapid  gun  fire  can 
be  on  the  order  of  6000  lbs,  to  8000  lbs.  Finally,  the  nearly 
uniform  undamaged  wing  drag  and  any  induced  drag  due  to  aileron 
control  deflection  needed  to  trim  the  aircraft  after  damage  will 
also  be  degrading  the  divergence  characteristics  of  the  wing. 
Possible  combinations  of  these  above  factors  adding  to  the  damaged 
induced  drag  forces  can  give  rise  to  a  degradation  of  the  diver¬ 
gence  speed  of  the  statistical  fighter  wing  to  within  the  flight 
envelope. 

A  further  consideration  of  the  evaluation  of  damage-in¬ 
duced  drag  modeling  can  be  found  in  appendix  A. 


3.  DAMAGE  INFLUENCE  ON  FLUTTER  AND  DIVERGENCE  SPEEDS 
OF  A  ZERO  SWEEP  STATISTICAL  FIGHTER  WING 

A  flutter  and  divergence  analysis  to  investigate  the  damage 
influence  on  the  statistical  fighter  wing  is  conducted  at  sea  level 
conditions  for  a  range  of  aerodynamic  damage  expressed  in  terms  of 
damage-induced  drag  levels.  During  the  studies,  three  aerodynamic 
aerodynamic  damage  configurations  are  considered.  In  one  case, 
aerodynamic  damage  is  assumed  localized  at  the  80%  semi-span  of 
the  wing,  while  in  the  second  case  this  damage  is  moved  to  90%  semi¬ 
span  position.  Finally,  a  third  case  of  60%  semi-span  location  is 
considered  to  estimate  the  tradeoff  between  aerodynamic  and  structural 
damage  location. 

For  a  cantilever  thin  wing  with  concentrated  aerodynamic  drag 
force,  the  basic  mathematical  model  employs  the  following  classical 
thin  beam  coupled  bending  and  torsional  equations  similar  to  those 
employed  in  lateral  buckling  studies  of  thin  beams*  (Appendix  B): 


H(x)  in  these  equations  represent  the  unit  step  function  and  is 
employed  here  to  more  compactly  write  the  beam  equations  that  are 
valid  to  the  right  and  left  of  the  concentrated  load. 


12 


13 


4  (ei  4>  * "  4 + s  4 

ax  ax  at  at 

2 

+  ADH(S-x)[(?-x)  -  2  ||]  -  L(w,B;  x,t) 

Sx 

l/r ,  ag^  ,  T  32g  ,  c  a20) 

ax  (G  eff  ax]  a  at2  Sa  ^2 

~2 

+  ADH(C-x)[(S-x)  -d-|]  -  M  (w,B;  x,t)  =  0 
3x^  ea 


=  0 


where 


AD  =  Cp  A  q 


*  1 

(Cp  =  constant  determined  from  experiment;  A  =  j  total  entrance 

plus  exit  hole  area) 

A  Galerkin-type  solution  of  these  equations  for  various 
levels  of  damage-induced  drag  is  employed  using 


“  =  V 

6  =  Boe 


iwt 


iwt 


wi  th 


Jo(x)  = 
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s0(x)  = 

and  h.,&.  defined  as  complex  with  f.(x),  F.(x)  determined  from 
i  J  1  J 

normal  vibration  modes.  Three  bending  f^(x)  and  three  torsional 
modes  F.j(x)  are  used  during  these  studies  for  different  classes 
of  structural  damage.  These  modes  are  presented  in  Tables  5 
through  14. 

The  appropriate  inertial  and  time-dependent  aerodynamic 
force  terms  are  also  included  in  the  analysis  thus  allowing 
a  complete  dynamic,  but  steady  state  stability  study.  For  the 
preliminary  studies,  incompressible  strip  theory  is  employed 
for  the  moderately  high  aspect  ratio  wings.  Finally,  a  constant 
static  parasite  drag  term  is  considered  to  be  the  significant 
or  first  order  aerodynamic  force  caused  by  the  damage.  These 
drag  estimates  are  deduced  from  Table  3  taken  from  Reference  9. 

The  results  of  drag  divergence  and  flutter  analyses  are 
presented  in  Tables  15  through  17  and  in  Figures  7  and  8.  The 
influence  of  aerodynamic  damage  alone  is  illustrated  in  Figure  7. 
This  damage  is  in  the  form  of  drag  only  with  100%  structural  inte¬ 
grity  assumed.  A  hypothetical  wing  tip  location  of  aerodynamic 
damage  site  (100%  semi-span)  is  also  presented  in  Figure  7.  This 
result  is  based  on  extrapolation  of  the  data  computed  for  60%, 

80%  and  90%  semi-span  locations  for  the  aerodynamic  sites.  The 


level  of  drag  is  seen  to  play  an  important  role  in  the  transition 
from  a  classical  flutter  critical  instability  for  the  undamaged 
wing  to  a  critical  divergence  type  of  instability  for  a  significant¬ 
ly  damaged  wing  subjected  to  high  drag  levels.  Interestingly  en¬ 
ough,  the  actual  flutter  speed  increases  as  a  result  of  the  added 
drag  term.  This  is  consistent  with  other  findings  in  the  litera¬ 
ture  for  moderate  aspect  ratio  undamaged  wings  (3,  4).  In  the 
present  study,  the  divergence  becomes  critical  due  to  the  high 
level  of  localized  drag  caused  by  the  damage.  The  near  frequency 
coalescence,  characteristic  of  classical  bending  torsional  flutter, 
is  eliminated  or  delayed  by  the  tendency  of  the  first  mode  fre¬ 
quency  branch  to  approach  zero  frequency,  thereby  promoting  diver¬ 
gence  in  the  first  mode  branch  prior,  to  the  occurence  of  flutter 
in  the  torsion  branch. 

In  spite  of  the  significant  reduction  in  the  critical 
divergence  speeds  of  the  wing  due  to  drag,  aerodynamic  damage 
alone  is  not  sufficient  to  reduce  these  critical  speeds  to  within 
the  flight  envelope  of  our  generic  fighter  wing  for  reasonable 
damage-induced  drag  levels.  Consequently,  several  structural 
damage  configurations  are  superimposed  upon  these  aerodynamic 
damage  cases  to  further  degrade  the  wing's  critical  speeds. 

These  results  are  presented  in  Tables  15  through  17  and  in  Figure 
8.  As  indicated  earlier,  the  structural  damage  is  imposed  by  lo¬ 
cally  reducing  the  torsional  rigidity  (CJ)  and  bending  stiffness 
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(El).  Without  aerodynamic  drag  imposed,  the  flutter  speed  is  seen 
to  decrease  due  to  the  loss  of  structural  stiffness.  These  re¬ 
sults  are  similar  to  the  findings  of  Reference  1.  For  the  larger 
reductions  of  El  and  GJ  at  stations  3,4,5,  the  critical  flutter 
speed  falls  into  the  flight  envelope  at  low  drag  level.  On  the 
contrary,  for  the  damage  at  station  5,6  the  divergence  speed  is 
more  critical  than  the  flutter  speed  at  moderate  to  high  damage- 
induced  drag  levels.  These  results  show  that  larger  structural 
damage  near  the  wing  root  will  produce  a  critical  flutter  condi¬ 
tion  while  aerodynamic  drag  rises  near  the  wing  tip  are  more  cri¬ 
tical  to  divergence.  Further  investigations  between  the  aerody¬ 
namic  damage  and  structural  damage  site  trade-offs  are  probably 
in  order  to  further  clarify  these  features. 

From  Figure  8,  for  aerodynamic  damage  at  the  30°'  semi-span, 
it  is  evident  that  for  maneuvering  flight  one  combination  of  aero¬ 
dynamic  and  structural  damage  equivalent  to  5%  of  the  wing  area  (see 
Figures  5&G)  and  producing  a  70%  reduction  in  torsional  and  bending 
stiffness  at  stations  4,  5  and  6  would  reduce  the  divergence  speed  of 
the  fighter  to  approximately  600  mph  or  to  within  its  flight  envelope 
Another  damage  alternative  occurs  for  an  aerodynamic  damage  at  90% 
semi-span  and  a  resulting  70%  reduction  in  torsional  stiffness  at 
stations  5  and  6  with  only  50%  reduction  in  bending  stiffness.  In 
this  case,  a  43,000  drag  load,  caused  by  a  4 V-%  hole  under  a  5  ’g' 
pullout,  will  reduce  the  divergence  speed  to  within  flight 
envelope.  Other  tradeoffs  are  also  possible,  as  is  evident  from 
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Figure  8.  Fighter  aircraft  wings  of  the  category  of  our  statistical 
model  are  probably  designed  to  a  load  factor  of  at  least  8.00  to  8.67 
with  a  1.5  margin  on  strength  under  these  conditions.  In  a  5  ' g ’  pull¬ 
out  maneuver,  therefore,  our  fighter  should  have  adequate  strength  even 
though  its  El  and  GJ  have  been  reduced  by  60%  to  70%  in  the  outer  wing 
panels  (10). 


4.  DAMAGE  INFLUENCE  ON  THE  FLUTTER 
AND  DIVERGENCE  OF  SWEPT  WINGS 

The  influence  of  ballistic  damage  on  the  flutter  and 
divergence  characteristics  of  swept  wings  is  also  investigated 
in  the  present  program.  The  study  employs  the  elastic  axis 
stiffness  distribution  of  the  statistical  fighter  wing.  The 
equations  of  motion  referenced  to  the  elastic  axis  can  be  re¬ 
written  as 


where  A  is  the  angle  of  swept  back  of  the  elastic  axis  and  n 
is  chordwise  location  of  the  damage-induced  drag. 
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Due  to  the  large  range  of  wing  sweep  considered,  the 
interaction  between  bending  and  twisting  deformations  greatly 
affects  the  modes  and  frequencies  of  vibration.  An  uncoupled 
bending  and  torsional  mode  vibration  analysis  is  no  longer  valid 
for  generating  assumed  modes  to  be  employed  in  a  Galerkin  analysis. 
Instead  of  using  the  uncoupled  modes,  a  Galerkin- type  solution  of 
the  flutter  equations  is  written  as: 

w  =  Jh.(x)  g.j(t) 
r 

e  =  j  e^x)  p.(t) 

where  g^(t)  and  p..(t)  are  the  generalized  coordinates  and 
h.j(x)  and  G^(x)  are  introduced  as  the  coupled  bending  and 
torsional  vibration  modes.  During  the  swept  wing  studies,  aero¬ 
dynamic  damage  is  located  at  90%  of  the  semi-span  and  different 
structural  damage  cases  are  imposed  at  Stations  5  and  6.  The 
six  coupled  modes,  used  for  this  flutter  study,  are  given  in 
Tables  18  through  21. 

The  numerical  studies  in  this  section  are  limit  checked 
against  the  results  of  Section  3  by  considering  the  special  case 
of  zero  wing  sweep.  A  comparison  of  these  computed  results  are 
illustrated  in  Table  22.  Column  'A'  here  represents  the  flutter 
and  divergence  results  for  zero  wing  sweep  taken  from  Section  3 


and  Table  17.  The  results  are  obtained  using  uncoupled  vibration 
modes  in  a  Galerkin  analysis.  Column  1 B *  is  a  solution  to  the 
same  problem  taken  from  the  studies  of  the  present  section  employ¬ 
ing  coupled  vibration  modes  in  a  Galerkin-type  solution  of  the 
flutter  equations.  The  divergence  speeds  are  seen  to  agree  with¬ 
in  1%  while  differences  in  computed  flutter  speeds  are  less  than 
5%. 

Wing  sweep  backs  of  20°  and  40°  are  investigated  for  a 
90%  semi-span  damage  site  location  while  a  swept  forward  wing  of 
-20°  is  also  checked  as  a  reference  configuration.  The  result 
of  the  study  for  aerodynamic  damage  only  with  n  =  0  is  illustrated 
in  Figure  9a  normalized  against  the  divergence  speed  of  the  res¬ 
pective  swept  wing  for  no  aerodynamic  damage.  A  review  of 
Figure  9a  suggests  that  a  critical  wing  sweep  exists  in  the  neigh¬ 
borhood  of  20°  that  is  most  susceptible  to  aerodynamic  damage. 

The  influence  of  two  chordwise  locations  of  the  aerodynamic  damage 
site  is  considered  for  a  sweep  of  +20°.  One  site  is  chosen 
forward  of  the  elastic  axis  0.833  ft  (g  <  0)  and  one  aft  of  the 
elastic  axis  by  0.833  ft  (n  >  0).  This  change  of  chordwise  damage 
site  location  is  found  to  produce  less  than  a  2%  difference  in 
the  flutter  and  divergence  speeds  between  the  two  cases.  Aero¬ 
dynamic  damage  induced  forward  of  the  clastic  axis  appears  to  be 
slightly  more  destabilizing.  The  influence  of  a  combined  aero¬ 
dynamic  and  structural  damage  investigation  is  also  illustrated  in 
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Figures  9b,  c  &  d.  The  results  demonstrate  that  structural  damage  de¬ 
grades  the  flutter  and  divergence  characteristics  of  swept-back  wings  to 
a  much  greater  extent  than  for  unswept  wings.  The  larger  the  sweep-back, 
the  larger  the  reduction  in  divergence  speeds  that  occur.  The  very  low 
divergence  speeds  associated  with  the  swept-forward  wing  are  also  il¬ 
lustrated  here.  The  apparent  lesser  influence  of  the  structural  damage 
on  swept-forward  wings  is  most  likely  due  to  the  already  inherently 
low  divergence  speeds  associated  with  these  configurations.  The  results 
of  the  trade-off  study  in  spanwise  damage  site  location  indicates  flutter 
again  may  be  the  crucial  instability  for  ballistic  damage  near  the  wing 
root,  while  drag  divergence  becomes  the  crucial  instability  for  damage 
near  the  wing  tip. 


5.  INFLUENCE  OF  DISTRIBUTED  AND  MULTIPLE  DAMAGE 
SITE  DRAG  EFFECTS  ON  FLUTTER  AND 
DIVERGENCE  OF  STRAIGHT  AND  SWEPT  WINGS 


The  investigations  of  Sections  3  and  4  can  be  generalized  to 
the  distributed  and  multiple  concentrated  drag  effects  illustrated 
in  Figure  10.  It  is  to  be  emphasized,  however,  that  the  actual  drag 
levels  and  spanwise  distribution  must  be  determined  from  wind  tunnel 
tests  as  no  theory  will,  at  present,  predict  these.  Undamaged  drag 
effects  can  also  be  included  in  this  generalization.  It  is  based 
upon  the  assumption  that  the  distributed  loading  can  be  "beamed"  or 
lumped  to  the  nodal  points  of  the  structure  as  is  done  in  any  standard 
finite  element  or  collocation  analysis.  In  the  case  of  a  series  of 
finite  concentrated  damage  sites  the  method  becomes  exact.  The  appro¬ 
priate  equations  are  readily  identified  as: 
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where  j  is  an  integer  and  N  is  the  number  of  concentrated  drag 
sites  or  nodal  points  used  in  discritizing  the  destributed  drag  loading. 
The  other  parameters  are  identified  in  Section  4. 


23 


To  illustrate  the  influence  of  distributed  drag  modeling  on 
wing  flutter  and  divergence  speeds  a  multiple  or  distributed  damage 
case  was  investigated  numerically.  The  results  of  this  study  are 
illustrated  in  Figure  11  for  four  cases  of  wing  sweep.  The  centroid 
of  the  loading  system  is  also  identified  in  these  figures.  Drag  levels 
at  Stations  5  and  6  were  taken  to  be  double  those  at  Station  4  for  this 
damage  configuration.  Structural  damage,  when  imposed,  was  considered 
occuring  only  at  Stations  5  and  6.  By  comparing  the  results  of  Figure 
11a  (zero  sweep  case)  with  the  data  of  Figure  8  for  damage  at  80% 
semi-span,  it  is  evident  that  the  general  trends  of  flutter  and  diver¬ 
gence  boundaries  with  drag  level  are  similar  for  distributed  and  con¬ 
centrated  drag  loadings.  A  further  review  of  these  figures  for  no 
structural  damage  indicates,  however,  that  replacing  the  distributed 
load  by  its  geometric  equivalent  concentrated  load,  acting  at  the 
centroid  of  the  distributed  drag  force,  predicts  too  low  a  flutter  and 
divergence  speed  estimates  for  this  case.  The  estimate  is  approximately 
20%  lower  than  it  should  be  for  divergence  speed  predictions.  This 
is  most  likely  due  to  the  fait  that  the  drag  is  actually  distributed 
over  1/3  of  the  semi -span  and  does  not  closely  approximate  a  concentrated 
loading  such  as  for  the  local i/ed  thru  hole  case.  Reviewing  the  other 
sweep  cases  indicates  as  before  that  structural  damage  has  a  more 
pronounced  influence  on  swept  back  configurations  than  on  swept  forward 
or  on  wings  with  zero  sweep. 


It  was  again  found  for  the  distributed  drag  cases  that  as  the 
damaged  induced  drag  level  increased,  the  flutter  speed  also  increased 
for  all  levels  of  imposed  structural  damage.  This  was  also  found  to 
be  the  case  for  the  concentrated  drag  configurations.  In  addition, 
wing  divergence  speeds,  as  in  the  concentrated  drag  case,  always  de¬ 
creased  with  increasing  damage  induced  drag  levels. 

In  many  failure  cases  an  aeroelastic  collapse  of  a  damaged 
wing  may  result  as  the  action  of  several  drag  factors.  As  an  illus¬ 
tration,  in  the  above  example  AD^  could  arise  from  external  store 
drag  or  from  the  recoil  of  a  rapid  firing  cannon  mounted  in  the  wing. 
The  damaged  induced  wing  drag  at  the  outer  tip  might  then  be  repre¬ 
sented  by  AD^  and  ADg.  The  recoil  from  a  20  min.  rapid  firing  cannon 
mounted  on  the  wing  might  be  on  the  order  of  4,000  to  8,000  pounds. 

For  the  above  simple  illustration  this  would  produce  total  chordwise 
loadings  on  the  wing  on  the  order  of  20,000  to  40,000  pounds.  A 
review  of  figures  5  and  6  for  sea  level  600  mph  flight  condition  indi¬ 
cate  that  such  damage  induced  drag  levels  of  16,000  and  32,000  pounds 
respectively  can  be  attributed  to  low  load  factor  flight  conditions. 

It  is  evident  from  this  study  that  distributed  drag  effects  allow  for  a 
combined  interaction  of  several  factors  which  may  require  an  extensive 
parameter  tradeoff  or  optimization  study  to  pin  down  the  extreme  cases 
of  interest.  The  above  illustration  merely  demonstrates  how  several 
interacting  drag  factors  can  be  integrated  into  the  analysis.  Drag 
divergence,  even  for  the  case  of  distributed  drag  influence,  is  still 
considered  to  be  a  critical  failure  mechanism  only  for  high  load  fac¬ 
tor  maneuvering  flight  conditions.  Further  parameter  tradeoff  studies, 
however,  are  still  needed  here. 


6.  COMPARISON  OF  THE  PREVIOUS  STRIP  THEORY  DRAG 
DIVERGENCE  MODELLING  WITH  A  LIFTING 
SURFACE  THEORY  MODELLING 
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The  strip  theory  aerodynamic  modelling  employed  for  lifting 
surfaces  investigated  in  the  present  and  earlier  studies  is  a  two- 
dimensional  approximation  to  a  three  dimensional  surface  configura¬ 
tion  valid  for  high  aspect  ratios.  In  view  of  this,  an  extension 
of  the  previous  drag  divergence  studies  to  lifting  surface  theory 
modelling  seems  necessary  to  further  evaluate  the  drag  divergence 
failure  mechanism  on  the  more  moderate  aspect  ratio  fighter  wings. 

Drag  divergence  investigations  were  conducted  employing  the 
kernel  function  lifting  surface  theory  of  Cunningham  outlined  in 
Reference  19.  These  studies  were  conducted  on  the  statistical  fighter 
wing  defined  in  Table  2,  and  on  a  modified  cantilever  model  of  the 
A-10  wing  illustrated  in  Figures  12  and  13.  The  results  of  the  study 
are  presented  below  and  plotted  in  Figure  9.  The  imposed  aerodynamic 
damage  induced  drag  is  located  at  the  90%  semi-span  position.  Structural 
damage  was  not  applied. 

Statistical  Fighter  Wing 

Damage  Indue ed  D rag_ 

Level  ~  lbs 

0 

25,919 
28,021 


Divergence Speed  mph 
Strip  Theory  Lifting  Surface 


1450  mph 


915 


1334  mph 
880 


25 


26 


Modified  A-10  Wing 


Damaged  Induced  Drag 
Level  ~  lbs 

0 

24,467 

25,509 


Divergence  Speed  mph 
Strip  Theory  Lifting  Surface 


890 

873 


920 

855 


In  all  of  the  above  cases  the  predicted  difference  in  divergence  speeds 
between  the  two  analyses  is  seen  to  be  less  than  9%  and  in  fact  less 
than  5"  in  nearly  all  cases  studied.  The  trends  in  the  statistical 
fighter  wing  furthermore  correclty  predict  the  changes  in  flutter  and 
divergence  properties  with  damage  induced  drag  that  were  calculated  for 
the  A-10  wing.  This  should  not  be  surprising  since  the  A-10  wing  is 
similar  in  design  and  structure  to  the  statistical  fighter  wing. 

In  summay  these  preliminary  studies  indicate  that  incompressible 
strip  theory  and  the  lifting  surface  theory  of  Reference  19  are  in  good 
agreement  when  it  comes  to  predicting  drag  divergence  failure  of  the 
statistical  fighter  wing  and  a  modified  A-10  wing. 
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7.  SIGNIFICANCE  OF  DAMAGED  INDUCED  UNSTEADY 
DRAG  EFFECTS  ON  PROMOTING  FAILURE  OF 
LIFTING  SURFACES 

It  has  been  shown  experimentally  in  wind  tunnel  tests  that  cer¬ 
tain  classes  of  damage  to  a  lifting  surface  can  produce  high  localized 
mean  drag  levels.  Such  high  drag  levels  occur,  as  shown  in  the  pre¬ 
vious  sections,  under  large  load  factor  maneuvering  flight  conditions 
and  can  lead  to  an  aeroelastic  failure  phenomenon  known  as  drag  diver¬ 
gence.  A  different  aeroelastic  failure  mechanism,  however,  can  also 
be  envisioned  even  for  non-maneuvering  flight  conditions  and  much  lower 
drag  levels  when  a  periodic  time  dependent  fluctuation  occurs  in  the 
damage  induced  drag  force.  (Figure  14).  Such  loadings  can  induce  both 
parametric  and  combination  resonances  on  a  wing  structure  (References 
14,  15,  16,  17,  18). 

Parametric  instabilities  or  resonances  of  structures  are  well 
studied  in  the  literature  and  occur  over  specific  ranges  of  frequencies 
of  forcing  function  (fluctuating  drag)  in  the  vicinity  of 

a)  »  2ton/k  k  =  1,2,3 

where  the  a>n's  are  the  characteristic  frequencies  of  the  structural 
system  (n  =  1,2,...).  For  non-conservative  systems  there  is  generally 
a  value  of  pulsating  forcing  level  below  which  parametric  resonance  is 
impossible,  irrespective  of  to  (References  14,  18).  So  called  combination 
resonance,  on  the  other  hand,  occurs  in  the  neighborhood  of 
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co  *  (u>n  +_  Wj)  k  =  1,2,... 

with  j  f  n.  Numerical  Floquet  analysis  has  been  shown  to  be  effective 
in  predicting  both  parametric  and  combination  resonances  (Reference  14). 

Figure  14  illustrates  the  implications  of  such  phenomenon  for 
amplitude  levels  of  time  dependent  drag  loadings  AD-j  and  mean  static 
loadings  ADq  that  may  be  much  less  than  the  critical  drag  divergence 
load  ADcr.  Unstable  regions  exist,  shaded  areas,  where  the  wing  under¬ 
goes  dynamic  oscillations  of  increasing  amplitude  produced  by  fluctuating 
drag  loads  that  are  appropriately  tuned  to  the  wing  natural  frequencies. 

The  appropriate  condition  of  tuning  is  roughly  given  by  the  above  expressions 
for  a).  It  is  not  difficult  to  envision  aerodynamic  damage  configurations 
that  give  rise  to  periodic  forces.  For  example,  the  blunted  trailing  edge 
of  a  wing  can  produce  periodic  shedding  and  air  loads  while  an  unstable 
shock  configuration  induced  by  damage  will  also  produce  oscillatory 
forces.  It  should  be  emphasized,  however,  that  an  adequate  definition  of 
such  time  dependent  forces  can  only  be  obtained  from  experimental  wind 
tunnel  studies  on  appropriately  damaged  models. 

Generally  the  nonconservative  nature  of  the  structural  damping 
forces  and  the  distributed  aerodynamic  forces  on  the  wing  will  preclude 
the  shaded  regions  from  touching  the  frequency  axis.  These  unstable 

regions  will  terminate  at  some  small  but  finite  AD./D  . x 

3  1  cr 

The  system  of  equations  developed  in  section  5  also  govern  the 
time  dependent  damaged  induced  drag  problem  when  the  drag  terms  AD.  are 

J 

assumed  to  be  time  dependent.  Since  a  significant  mean  drag  term  ADqj 
may  develop  due  to  damage  in  addition  to  a  fluctuating  drag  component, 
the  drag  AD^  is  expressed  as 
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ADj  =  AD0j  +  AD,J  F(t)  , 

F(t)  =  F(t  +  T) 

T  =  Period  of  Motion 

where  the  function  F(t)  as  determined  from  wind  tunnel  tests,  can  be 
identified  in  terms  of  its  Fourier  components..  The  solution  of  the 
resulting  partial  differential  equations  with  variable  but  periodic 
coefficients  is  then  reduced  to  a  system  of  Hill  equations  thru  a  Galerkin 
approximation  of  the  spacial  variables  as  outlined  in  Section  3. 

Numerical  Floquet  analysis  can  then  be  effectively  employed  to  solve  this 
system  of  Hill  equations  to  determine  the  regions  of  both  parametric  and 
combination  resonances  illustrated  by  the  shaded  regions  of  Figure  14. 

To  illustrate  the  occurance  of  this  phenomenon,  preliminary 
numerical  solutions  were  obtained  for  the  statistical  fighter  wing 
studied  earlier  in  this  report.  For  a  drag  force  fluctuating  harmonically 
with  frequency  a  two  mode  Galerkin  analysis  was  employed  utilizing 
the  first  bending  and  torsion  modes  of  the  straight  elastic  axis  fighter 
wing.  The  resulting  coupled  Hill's  equations  were  then  solved  for  small 
values  of  ADQ/ADcr  employing  the  perturbation  procedure  of  Reference  17. 
Since  Hsu's  scheme  is  restricted  to  only  a  special  class  of  nonconservative 
systems  the  distributed  aerodynamic  lift  and  moment  where  initially  ex¬ 
cluded  and  only  the  time  dependent  drag  terms  included  in  the  analysis 
as  the  first  order  fluid  effect.  The  results  of  the  study  are  illustrated 
in  Figure  15  where  the  shaded  areas  indicate  the  regions  of  instability. 
Should  the  periodic  forcing  frequency  of  the  damaged  induced  drag  force 
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fall  within  the  shaded  regions  of  the  plot  the  wing  structure  will  become 
dynamically  unstable.  That  is,  its  oscillatory  motion  will  increase 
exponentionally  with  time. 

The  study  for  the  inclusion  of  distributed  lift  and  moment  as 
well  as  for  the  fluctuating  damage  induced  drag  loading  must  include  a 
transient  aerodynamic  force  modelling.  This  has  been  initially  approxi¬ 
mated  by  a  quasi-steady  strip  theory  for  subsonic  Mach  numbers  and  a 
piston  theory  approximation  for  supersonic  Mach  numbers.  A  classical 
Galerkin  solution  is  then  employed  along  with  an  eighth  order  variable 
step  Runge-Kutta  integrator  in  a  numerical  Floquet  analysis  to  examine 
the  stability  of  the  damaged  wing  configuration.  These  results  are 
currently  under  development  but  preliminary  studies  here  are  in  basic 
agreement  with  Figure  15  except  near  ADQ/ADcr  0  where  aerodynamic 
damping  effects  eliminate  any  possible  unstable  regions. 


8.  CONCLUSIONS 


Aeroelastic  failure  modes  due  to  warhead  damage  have  been  investi¬ 
gated  for  several  lifting  surface  configurations.  The  damage  model 
includes  both  structural  and  aerodynamic  modifications  to  the  lifting 
surface.  Structural  damage  is  reflected  as  reductions  of  the  stiffness 
distributions  El  and  GJ  in  the  region  of  the  damage  while  aerodynamic 
damage  is  included  to  first  order  in  the  form  of  a  drag  rise  localized 
at  the  damage  site.  Classes  of  damage  producing  relatively  small  steady 
state  drag  rises  but  significant  reductions  in  El  and  GJ  inboard 
on  the  wing  degrade  the  flutter  speed  of  a  statistical  fighter  wing  to 
within  its  flight  envelope.  On  the  other  hand,  for  structural  damage 
in  the  outer  wing  panels  which  produces  high  steady  state  drag  rise, 
divergence  becomes  the  critical  mode  of  instability.  Under  high  load 
factor  maneuvers,  this  divergence  speed  can  drop  to  within  the  flight 
envelope  of  the  aircraft. 

The  influence  of  warhead  damage  on  the  aeroelastic  response  of 
swept  wings  is  carried  out  over  a  range  of  swept  back  and  swept  for¬ 
ward  planforms.  The  study  indicates  that  aerodynamic  damage  in  the  form 
of  drag  rise  is  the  most  critical  to  wings  with  approximately  a  20° 
swept-back  planform.  Structural  damage  also  appears  more  critical  to 
swept  back  wing  configurations.  The  greater  the  wing  sweep-back,  the 
greater  the  influence  of  structural  damage  in  reducing  the  flutter  and 
divergence  speeds  of  the  wing.  Additional  trade-off  studies  for  damage 
site  locations  on  swcpt-back  wings  are  desirable. 
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The  results  of  the  present  study  are  consistent  with  the  findings 
of  Reference  1  when  aerodynamic  damage  is  ignored  and  only  structural 
damage  is  considered. 

An  evaluation  of  the  drag  divergence  modeling  employing  an  incom¬ 
pressible  strip  theory  aerodynamic  approximation  provided  drag  divergence 
boundaries  within  5 %  of  those  predicted  by  a  more  refined  lifting  sur¬ 
face  theory.  In  the  case  of  the  statistical  fighter  wing  the  boundaries 
predicted  by  strip  theory  were  found  to  be  3  to  91  higher  in  velocity 
than  those  predicted  by  lifting  surface  theory;  while  for  a  modified 
A-10  wing  the  divergence  speeds  were  approximately  3  to  5%  below  those 
estimated  by  the  lifting  surface  aerodynamic  theory.  In  the  summary, 
the  aeroelastic  drag  divergence  modeling  employing  incompressible  strip 
theory  is  providing  an  accurate  estimate  of  the  drag  divergence  features 
for  a  damaged  wing. 

A  second  aeroelastic  failure  mechanism  can  be  identified  for 
damaged  wings  when  a  periodic  chordwise  force  is  generated  by  the  aero¬ 
dynamic  damage.  Under  these  conditions  "so  called"  parametric  and 
combination  resonance  or  instability  has  been  shown  to  occur  when  the 
periodic  chordwise  force  is  tuned  to  the  structural  resonant  frequencies 
of  the  wing.  The  magnitude  of  the  drag  forces  producing  this  instability 
can  be  much  less  than  those  causing  drag  divergence. 

Since  the  validity  of  the  proposed  damage-induced  aeroelastic 
failure  mechanism  are  strongly  dependent  upon  an  accurate  assessment 
of  aerodynamic  damage  in  the  form  of  damage-induced  drag  increments, 
more  experimental  wind  tunnel  studies  over  the  results  of  Reference  11 
appear  necessary.  Finally,  further  refinements  in  the  aeroelastic 
modeling  are  also  recommended  to  allow  for  the  time  dependent  character 
of  the  warhead  blast  loads. 
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a)  Finite  element  lumped  mass  modeling  of  wing  structure  along 
its  elastic  axis 


NOTE:  Locations 
of  the  two  aero¬ 
dynamic  damage 
sites  investigated 
CD-  80%  semi -span 
0-  90%  semi -span 


refere  nee 
axis 


b)  Two  spar  stressed  skin  semimonocoque  wing  structure 


trailing 

edge 


c)  Strip  theory  aerodynamic  modeling 

-b  0  ba  be  be  +b  x=bx 


AEROELASTIC  MODELING  EMPLOYING  AN  ELASTIC 
AXIS  BEAM  TYPE  STRUCTURE  AND  A  STI  IP 
THEORY  AERODYNAMIC  APPROXIMATION 


Figure  1 


I 


IDEALIZED  AERODYNAfl] 


Fi gure 


YP3  00  at  wing  root 


YP3 1.00  at  wing  tip 


YP=,71  .79 


YP  3  .57 

a  =  Undamaged 
•  3  Damaged 


.28 

•3'64  ^f1 68.96 

ACP 

Piot  Locations 

YP  3  .74 

a  =  Undamaged 
•  =  Damaged 

.17  33  50  67  83  100 

x/c 


.33  .50  87  S3  l.OC 

X/C 


YP  3  .38 


a  =  U'ndomoged 
•  Damoged 


.17  33  50  67 

X/C 


THEORETICALLY  PREDICTED  PRESSURE  DISTURBANCES  PRODUCED  BY  A  THRU-HOLE 
TYPE  DAMAGE  (SHADED  REGION)  ON  A  MODERATE  ASPECT 
RATIO  LIFTING  SURFACE  WITH  TAPER 

Figure  3 


ATIVE  AND  UPPER 
FOR  VARIOU: 


AD~Damage  Induced  Drag  ~  LB  x  10 
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j 


^  Concentrated  Drag  100% 
^  Semi -span  -  Sea  Level 


True  Speed  Up~Critical  Divergence  Speed~mph 
Up~Critical  Rutter  Speed  ~mph 


INFLUENCE  OF  SPANWISE  LOCATION  OF  AERODYNAMIC  DAMAGE 
SHE  (NO  STRUCrURAL  DAMAGE)  ON  DRAG 
DIVERGENCE  AND  FLUTTER 


Figure  7 


AD~Damage  Induced  Drag  ~  LB  x  I0-3 
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Aerodynamic  Damage  at  60% 
Semi- Span /Structural  Damage 


0  400  800  1800  1600 

True  Speed  UD' 

Up’ 


Aerodynamic  Damage  at  80% 
Semi-Span/Structural  Damage 
at  Stations  4, 5  and  6 


Aerodynamic  Damage  at  90% 
Semi-Span /Structural  Damage 
at  Stations  5  and  6 


Aerodynamic  Damage  at  95% 
Semi -Span /Structural  Damage 
at  Stations  5  and  6 


Decrease 


El 

30 

40 

50 

70 


No  Structural 
Damage  Only 
Aerodynamic 
Damage 

—  U0  Divergent  Speed 
• . Up  Flutter  Speed 


0  400  800  1200  1600  1800 

■Critical  Divergence  Speed~mph 
■Critical  Flutter  Speed  ~mph 


INFLUENCE  OF  SPANWISE  I.OCAIION  OF 
AERODYNAMIC  AMD  STRUCTURAL  DAMAGE  SITE 
ON  DRAG  DIVERGENCE  AND  FLUTTER 


Figure  8 


AD~Damage  Induced  Drag  ~ LB  x  10" 
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A- 10  DRAG  DIVERGENCE  MODES 
Figure  12 


Mode  6 


AD/ADcr 


INFLUENCE  OF  TIME  DEPENDENT  DAMAGE  INDUCED  DRAG 
FORCE  ON  PROMOTING  PARAMETRIC 
TYPE  FUJ  I  TER  INSTABILITIES 


Figure  14 


PARAMETRIC  INSTABILITIES  INDUCED  IN  STATISTICAL  FIGHTER 
BY  DAMAGE-INDUCED  DRAG,  PERIODIC  IN  TIME 


TABLE  1 


STATISTICAL  FIGHTER  PROPERTIES  DETERMINED  AS 
GEOMETRIC  MEANS  OF  THE  FOURTEEN  DIFFERENT 
FIGHTER  AIRCRAFT 


ce  Designation 

Manufacturer 

P-51  A 

North  American 

XP-60 

Curtiss-Wright 

XP-84 

Republic 

XFJ-1 

North  American 

XF8F-1 

Grumman 

XF9F-2 

Grumman 

F9F-3 

Grumman 

XFD-1 

McDonnel 1 

F2H-1 

McDonnell 

XF2D-1 

McDonnell 

XF6F-1 

Grumman 

XBT2D-1 

Douglas 

XF3D-1 

Douglas 

XF6U-1 

Chance-Vought 

Parameter  Xg  of  typical  fighter  determined  as  geometric  mean  of 
similar  parameters  X^,  X.,  ,  etc.  of  above  aircraft.  That  is, 

x3  =  (X,.  X2  •  X3  ....xra),/ra 

where  X^ ,  X^,  - Xm  represent  the  set  of  variates.  The  speci¬ 

fic  typical  wing  parameters  determined  were 

Height  •  Chord  length 

Static  moment  about  a  reference  axis  •  Chordwise  center  of  gravity 
Pending  area  moment  of  inertia  position 

Mass  moment  of  inertia  about  a  ref-  •  Chordwise  elastic  axis  posi- 
erence  axis  tion 

Co- efficient  of  torsional  rigidity 


CO  VO  VO 


CO 

o 


Mcoinroioojcoin 


ro  vo  m  o 


lO  CO  vO  CO  O  *3-  LT) 


CO  VO  VO 


co  vo  CO 


AIRPLANE  473  5.26  352  3.91  247  2.74  138.  1.53 


TABLE  II  Cont'd. 


Semispan  -  235  in. 


TABLE  3. 

KEY  TO  GENERALIZED  DAMAGE  CLASSES 
(Taken  from  Reference  9) 
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Dashed  lines  in  the  sketches  below  indicate  regions  of  possible  hole  damage. 
Projectile  penetrates  both  regions  of  each  sketch.  These  dashed  lines  do  rot  indi¬ 
cate  total  skin  damage,  but  locations  of  a  family  of  possible  projectile  holes. 

Poles  of  class  [  S  F  damages  have  projecting  skin;  all  other  damages  are  flush  holes. 


Damage 

Classes 


Wing  Section 


C,D 


e,g,h  <;; 


Damage  Classes 

Hole  Sire* 

A. C.E.F.G.J 

B, D,  H,  K 

,15c 

10c 

Diameter  of  an  equal  area  circular  hole 


Jligest  of  Drag  Increases  Due  to  Damage 


Damage 

Class 


'  Dial  .ter  of  on  equal  area  circular  hole 


i 
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TABLE  4 

DAMAGE  INDUCED  DRAG  LEVELS  FOR  VARIOUS 
MANEUVERING  FLIGHT  CONDITIONS 


LOAD  FACTOR 

CL 

★ 

CD 

AD  LBS 

n 

upper  & 
lower  bounds 

2%  hole 

3%  hole 

5%  hole 

1 

0.143 

1 .10 

2.60 

3,242# 

to 

7,660# 

4,863# 

to 

11,492# 

8,150# 

to 

19,150# 

5 

0.713 

3.54 

7.48 

10,433# 

to 

22,040# 

15,649# 

to 

33,062# 

26,080# 

to 

55,100# 

6 

0.856 

4.70 

9.73 

13,851# 

to 

28,667# 

20,777# 

to 

43,000# 

34,628# 

to 

71  ,667# 

6.5 

0.927 

5.30 

11.0 

15,619# 

to 

32,400# 

23,429# 

to 

48,620# 

39,048# 

to 

81 ,033# 

600  mph  true  maneuvering  speed 
20,000  ft.  combat  altitude 
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TABLE  5 


FIRST  SIX 

UNCOUPLED  VIBRATIONAL 

MODES 

NO 

STRUCTURAL  DAMAGE 

I . 

Bending  Modes 

a),  =  55.849 

(jj  2  =  188.314 

“3  = 

1 . 

1  .  7857E-2 

-8. 4490E-2 

2 

2  . 

8. 1406E-2 

-3.111 3E-1 

6 

3. 

2 .0923E-1 

-5 . 4403E-1 

3 

4. 

4. 1  146E-1 

-5.1 329E-1 

-7 

5. 

6.8315E-1 

1 . 4467E-2 

-9 

6. 

1  .OOOOE  +  O 

1 .OOOOE+O 

1 

II 

Torsional  Modes 

u  i  =  188.613 

u2  =  342.126 

w3  = 

1  . 

9.6133E-2 

-1 .8214E-1 

2 

2. 

2 . 2696E-1 

-3. 3226E-1 

2 

3. 

4.  1  644E-1 

-3. 2308E-1 

-2 

4. 

6 . 4  0  3  2  E  -  1 

-1 . 6795E-2 

-5 

5. 

8.  395  1  E-1 

4. 7194E-1 

-1 

6. 

1  .OOOOE  +  O 

1 .OOOOE+O 

1 

450.21  7 
( rad/ sec ) 

4423E-1 

1549E-1 

0448E-1 

5282E-1 

5454E-1 

OOOOE+O 


504.685 
( rad/sec) 

41 34E-1 

2141E-1 

.  4823E-1 

9986E-1 

4908E-1 

0000E+0 
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TABLE  6 

FIRST  SIX  UNCOUPLED  VIBRATIONAL  MODES 
El- 30% ,  G J -50%  DECREASED  0  STA.  3,4,5 


I . 

Bending  Modes 

(x) 

=  52.640 

ix>2  =  170.078 

“3 

1 . 

1  .  5054E-2 

-7 . 9392E-2 

-2 

2. 

6.8823E-2 

-3.0221 E - 1 

-3 

3. 

1  .8015E-1 

-5.6571E-1 

-3 

4. 

3. 7690E-1 

-5 . 8070E-1 

6 

5. 

6.6068E-1 

-4 . 6049E-2 

1 

6. 

1  .OOOOE  +  O 

1  .OOOOE  +  O 

-9 

II. 

Torsional  Modes 

Fj(x) 

uj  =  157.419 

a>2  =  292.544 

“c 

1  . 

4.8591E-2 

-  1  . 7891  E-1 

1 

2. 

1  .  181 7E-1 

-3.6344E-1 

2 

3. 

2 . 6 5 72 E -  1 

-5  .  1  796E-1 

-1 

4. 

5 . 6092 E  -  1 

-  2  .  5  3 1  2  E  -  1 

-7 

5  . 

8. 3818E-1 

4.41  14E-1 

-1 

6. 

1  .OOOOE  +  O 

1 .OOOOE+O 

1 

407.619 

(rad/sec) 

. 0984E - 1 

. 5680E-1 

.  8472E-1 

.  87  30E- 1 

.  OOOOE  +  O 

.5720E-1 


421 . 626 
(rad/sec) 

. 991 7E-1 

. 8  36  7  E -  1 

. 2760E-2 

. 7386E-1 

.  6086E- 1 

.  OOOOE  +  O 


TABLE  7 


FIRST  SIX 

UNCOUPLED  VIBRATIONAL 

MODES 

£1-40%, 

GJ -60%  DECREASED  G>  STA. 

3,4,5 

I . 

Bending  Modes 

u,  =  51.127 

gj  2  =  1  63.074 

“3  = 

1 . 

1 . 3863E-2 

-7. 7807E-2 

-1 

2  . 

6 . 3465E-2 

-2 . 9988E- 1 

-5 

3. 

1  . 6  742  E - 1 

-5 . 762 1 E - 1 

-4 

4. 

3.6138E-1 

-6. 1283E-1 

6 

5. 

6 . 5057E-1 

-7.5155E-2 

1 

6. 

1 .0000E+0 

1 .0000E+0 

-9 

I  I 

.  Torsional  Modes 

u)  -j  =  148.836 

=  286.953 

w3  = 

1  . 

4. 2014E-2 

-1 . 9  66 1  E -  1 

1 

2. 

1 . 0290E-1 

-3.9334E-1 

2 

3. 

2.4635E-1 

-5 . 9470E-1 

3 

4. 

5 . 8476E-1 

-2 . 7798E-1 

-8 

5. 

8.5998E-1 

4. 7952E-1 

-4 

6. 

1 . OOOOE+O 

1 .OOOOE+O 

1 

389.939 

(rad/sec) 

. 9324E-1 

. 3870E-1 

. 09 1 4E- 1 

. 4797E-1 

. 0000E+0 

. 1886E-1 

406.357 

(rad/sec) 

. 8453E-1 

. 781 7E-1 

. 0931 E-2 

. 3054E-1 

. 3754E-2 

. 0000E+0 
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TABLE  8 


FIRST  SIX 

UNCOUPLED  VIBRATIONAL 

MODES 

EI-50%, 

G J -  70%  DECREASED  0  STA. 

3,4,5 

I . 

Bending  Modes 

=  49.225 

=  155.417 

w,  =  369.508 
(rad/sec) 

1 . 

1  . 2479E-2 

-7 . 6342E-2 

-1 . 7589E-1 

2. 

5 . 7220E-2 

-2 . 982 1 E - 1 

-5.0725E-1 

3. 

1 . 5234E-  1 

-5 . 8956E- 1 

-4 . 3585E- 1 

4. 

3 . 4278E-I 

-6 . 5348E-1 

6 . 0486E- 1 

5. 

6 . 3849  E - 1 

-1 . 1 22 1 E - 1 

1 .OOOOE+O 

6. 

1 . OOOOE+O 

1 .OOOOE+O 

-8.851 7E-1 

1 1 

Torsional  Modes 

uj  =  132.217 

U2  =  264.400 

a>7  =  361 . 831 
°  (rad/sec) 

1  . 

2 . 8352E-2 

-1 . 5937E-1 

2 . 00  50  E - 1 

2  . 

7.0314E-2 

-  3 . 4  0  2  1  E  - 1 

3. 4743E-1 

3. 

1 . 8395E-1 

-6. 0422E-1 

2.4247E-1 

4. 

5 . 2435E-1 

-5 . 0652  E -  1 

-9.5912E-1 

5. 

8. 5386E-1 

4. 1558E-1 

-9  .  4496E-2 

6. 

1 .0000E+0 

1 . OOOOE+O 

1 .OOOOE+O 
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TABLE  9 


FIRST  SIX 

UNCOUPLED  VIBRATIONAL 

MODES 

m 

*—* 

i 

co 

O 

G J -50%  DECREASED  @  STA. 

4,5,6 

I . 

Bending  Modes 

i  =  54.586 

u2  =  173.882 

w,  =  416.531 

(rad/sec) 

1 . 

1  . 6462E-2 

-7 . 5298E-2 

-2. 1644E-1 

2  . 

7.5150E-2 

-2.841 1 E - 1 

-5 . 7309E-1 

3. 

1  .9354E-1 

-5 . 2369E- 1 

-3.6521E-1 

4. 

3.8446E-1 

-5. 5831E-1 

6. 3126E-1 

5. 

6.6153E-1 

-5 . 4452E-2 

1 . 0000E+0 

6. 

1  .  0000E  +  0 

1 .0000E+0 

-8 . 8978E-1 

1 1 

Torsional  Modes 

u1  =  169.266 

U)2  =  288.424 

oi,  =  426.189 

(rad/sec) 

1  . 

5.2174E-2 

- 1 . 4506E -  1 

1  .9416E-1 

2. 

1  .2556E-1 

-2 . 9697E-1 

2. 7157E-1 

3. 

2. 3897E-1 

-3.9974E-1 

4 . 2939E-2 

4. 

4. 4240E-1 

-3.0187E-1 

-5 . 9255E-1 

5. 

7.4149E-1 

2.4940E-1 

-6. 3887E-1 

6. 

1 . 0000E+0 

1 . 0000E+0 

1 .0000E+0 

TABLE  10 


FIRST  SIX 

UNCOUPLED  VIBRATIONAL 

MODES 

El -40% , 

GJ-60%  DECREASED  @  STA. 

4,5,6 

I . 

Bending  Modes 

u»1  =  53.959 

u>2  =  1  6  7.666 

=  403.325 
(rad/ sec) 

1 . 

1 . 5803E-2 

-7. 1950E-2 

-1  .  9  7  4  6  E  - 1 

2. 

7.2194E-2 

-2. 7435E-1 

-5. 3308E-1 

3. 

1 .8611E-1 

-5. 1696E-1 

-3. 7068E-1 

4. 

3.7132E-1 

-5 . 7898E-1 

5 . 5428E-1 

5. 

6 . 5072E-1 

-8.6361E-2 

1 . OOOOE+O 

6. 

I  . 0000E+0 

1  . OOOOE  +  O 

-8. 1275E-1 

II 

.  Torsional  Modes 

u1  =  160.891 

oj2=  274.021 

w,=  400.230 
(rad/sec) 

1  . 

4. 0782E-2 

-1 . 4069E-1 

1  .  7509E-1 

2. 

9 . 8884E-2 

-2  .  9552E-1 

2 . 6964E-1 

3. 

1  .  9089E-1 

-4.2316E-1 

1  .  2086E-1 

4. 

3.  7865E-1 

-3.9885E-1 

-4. 9533E-1 

5  . 

7 . 0805 E - 1 

1 . 531 3E-1 

-8 . 0663E-1 

6. 

1  .0000E+0 

1  .0000E+0 

1 . 0000E+0 

■MH* 


TABLE  11 


FIRST  SIX 

UNCOUPLED  VIBRATIONAL 

MODES 

EI-50%, 

GJ -70%  DECREASED  @  STA. 

4,5,6 

I . 

Bending  Modes 

a)  i  =  5  3.143 

to  2  =  160.410 

ui,  =  388 . 571 

0  (rad/sec) 

1 . 

1  . 4980E -2 

-6. 8511E-2 

-1 . 8581 E-l 

2. 

6 . 8498E-2 

-2 . 6448E-1 

-5. 1259E-1 

3. 

1 .7681  E-1 

-5.1102E-1 

-3. 9055E-1 

4. 

3.5461E— 1 

-6.0478E-1 

4.9183E-1 

5. 

6 . 3673E-1 

-1 .2608E-1 

1 .OOOOE+O 

6. 

1  .OOOOE+O 

1 .OOOOE+O 

-7 . 6730E-1 

II 

.  Torsional  Modes 

u1  =  148.750 

u2  =  257.867 

u,  =  365.769 

(rad/sec) 

1  . 

2 . 8938E-2 

-1 . 3573E-1 

-1 .5059E-1 

2. 

7 . 0883E-2 

-2 . 9280E-1 

-2 . 5832E-1 

3. 

1  .  3944E-1 

-4 . 4553E-1 

-2.01 29E-1 

4. 

3.0350E-1 

-5.2027E-1 

3.1 150E-1 

5. 

6 . 6726E- 1 

4. 4689E-5 

1 . OOOOE+O 

6. 

1  .OOOOE  +  O 

1 .OOOOE+O 

-9.91 49E-1 
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TABLE  12 


FIRST  SIX 

UNCOUPLED  VIBRATIONAL 

MODES 

El -30% , 

G J  -50%  DECREASED  @  STA. 

5,6 

I . 

Bending  Modes 

w-|  =  55.593 

u2  =  183.106 

w3 

1 . 

1  .7496E-2 

-7 . 8475E-2 

-2 

2. 

7 . 9785E-2 

-2.9138E-1 

-5 

3. 

2.0516E-1 

-5 . 1869E-1 

-3 

4. 

4 . 0470E- 1 

-5. 1456E-1 

6 

5. 

6.7315E-1 

-  3 . 6596E-2 

1 

6. 

1 .0000E+0 

1 . 0000E+0 

-9 

II. 

Torsional  Modes 

to  1  =  181.590 

u2  =  302.679 

u3  = 

1  . 

6.8795E-2 

-1 . 1984E-1 

2 

2. 

1  . 6359E-1 

-2 . 3870E-1 

3 

3. 

3.0439E-1 

-2 . 9968E-1 

-2 

4. 

4. 7818E-1 

-1 .9779E-1 

-6 

5. 

7 . 0248E-1 

1 . 7338E-1 

-7 

425.207 

(rad/sec) 

. 1294E-1 

. 5685E - 1 

. 3499E-1 

.  1  895E-1 

.  0000E+0 

. 1 954E-1 


443.677 

(rad/sec) 

.  4890E-1 

. 2  32 1  E - 1 

.  0751  E-2 

.0198E-1 

. 761 3E-1 


6. 


1 . 0000E+0 


1  . 0000E  +  0 


1 . 0000E+0 


TABLE  13 

FIRST  SIX  UNCOUPLED  VIBRATIONAL  MODES 
El -40%  ,  G J - 60%  DECREASED  0  STA.  5 ,6 
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I. 

Bending  Modes 

u»1  =  55.471 

w2  =  180.628 

w,  =  414.239 
(rad/ sec) 

1 . 

1  . 7324E-2 

-7.5778E-2 

-1 .9813E-1 

2. 

7.9015E-2 

-2 . 8248E-1 

-5.2652E-1 

3. 

2 . 0322E-1 

-5.0710E-1 

-3.4169E-1 

4. 

4 .  OOOOE-l 

-5. 1463E-1 

5.5501E— 1 

5. 

6.6819E-1 

-6 . 0867E-2 

1 .OOOOE+O 

6. 

1  -OOOOE+O 

1  .OOOOE+O 

-8.6903E-1 

II  . 

Torsional  Modes 

u>i  =  178.142 

«2  =  287.737 

u>3  =  426.199 

(rad/sec) 

1  . 

5 . 8882E-2 

-1 .0887E-1 

-2 . 5879E-1 

2. 

1  . 4050E- 1 

-2.2316E-1 

-3. 6197E-1 

3. 

2.6315E-1 

-3.01 34E-1 

-5.71 88E-2 

4. 

4.1 754E-1 

-2.4928E-1 

5. 6523E-1 

5. 

6.4208E-1 

6 . 6228E-2 

1  .OOOOE+O 

6. 


1 . OOOOE+O 


1 .OOOOE+O 


-9. 5358E-1 


TABLE  14 


FIRST  SIX 

UNCOUPLED  VIBRATIONAL 

MODES 

EI-50%, 

GJ - 70%  DECREASED  @  STA. 

5,6 

I . 

Bending  Modes 

to,  =  54.878 

a>2  =  1  68.958 

“3  = 

1 . 

1 . 6500E-2 

-6.4717E-2 

-1 

2. 

7.5313E-2 

-2 . 4579E-1 

-4 

3. 

1  .  9  3  9  2  E  - 1 

-4. 5888E-1 

-3 

4. 

3 . 822 1 E - 1 

-5. 1422E-1 

3 

5. 

6 . 4323E-1 

-1 . 7061 E-l 

1 

6. 

1  .OOOOE+O 

1 .OOOOE+O 

-6 

II 

Tors i onal  Modes 

=  172.316 

a)2  =  268 . 505 

“3 

1  . 

4. 5928E-2 

-1  .0244E-1 

-2 

2. 

1 . 1021E-1 

-2  .  1  720E-1 

-3 

3. 

2 . 0863E- 1 

-3.1783E-1 

-1 

4. 

3.3641E-1 

-3.2135E-1 

3 

5. 

5 . 5348E-1 

-8. 41 58E-2 

1 

6. 

1 .OOOOE+O 

1  .OOOOE  +  O 

-6 

372.134 

(rad/sec) 

.  5806E-1 

. 4640E-1 

.7213E-1 

. 4990E-1 

.  OOOOE+O 

. 88  70  E-  1 


406.1  36 
( rad/ sec) 

. 1 376E-1 

. 2250E-1 

. 2443E-1 

. 8566E- 1 

.OOOOE+O 

.  7547E-1 
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TABLE  15 

CRITICAL  SPEEDS  FOR  DAMAGED  WING 
AERODYNAMIC  DAMAGE  SITE  AT  60%  SEMI-SPAN 


No  Structural 

Damage 

di vergence 

D 

f  1  utter 

D 

speed 

lbs 

speed 

lbs 

mph 

mph 

1470 

0 

950 

0 

1239 

25,678 

968 

15,673 

1134 

43,021 

990 

32,789 

983 

80,735 

1133 

107,254 

870 

126,609 

— 

— 

EI-30%,  GJ - 50% 

Decreased  @  STA. 

3,4,5 

divergence 

D 

flutter 

D 

speed 

lbs 

speed 

1  bs 

mph 

mph 

1165 

0 

734 

0 

1025 

17,574 

748 

9,359 

837 

58,592 

820 

56,236 

748 

93,689 

— 

— 
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TABLE  15  CONT'D 


3.  El -40% ,  GJ -60%  Decreased  0  STA.  3,4,5 


d i vergence 
speed 
mph 

D 

lbs 

flutter 

speed 

mph 

D 

lbs 

1098 

0 

690 

0 

969 

15,715 

701 

8,224 

794 

52,700 

768 

49,305 

711 

84,464 

— 

— 

EI-50%,  G J - 70% 

Decreased  @  STA. 

3,4,5 

divergence 

D 

flutter 

D 

speed 

lbs 

speed 

lbs 

mph 

mph 

957 

0 

588 

0 

865 

12,515 

596 

5,941 

723 

43,682 

640 

34,228 

651 

70,912 

— 

— 
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TABLE  16  .  ; 

CRITICAL  SPEEDS  FOR  DAMAGED  WING  j 

AERODYNAMIC  DAMAGE  SITE  AT  80%  SEMI-SPAN  f 

\ 

] 

1.  No  Structural  Damage  j 


di vergence 
speed 
mph 

D 

1  bs 

flutter 

speed 

mph 

D 

lbs 

1470 

0 

950 

0 

1151 

22,160 

980 

16,065 

1038 

36,050 

1040 

35,91 1 

887 

65,830 

— 

— 

785 

102,950 

— 

— 

EI-30%,  GJ -50% 

Decreased  @  STA. 

4,5,6 

di vergence 

D 

flutter 

D 

speed 

lbs 

speed 

1  bs 

mph 

mph 

1234 

0 

784 

0 

1010 

17  ,060 

808 

10,918 

794 

52,730 

— 

— 

708 

83,850 

— 

— 

643 

138,320 

— 

— 

5 
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TABLE  16  CONT'D 

5.  E 1-70% ,  GJ - 70%  Decreased  @  STA.  4,5,6 


di vergen  ce 
speed 
mph 

D 

lbs 

f 1  utter 
speed 
mph 

D 

lbs 

1040 

0 

649 

0 

865 

12,256 

668 

7,309 

687 

39,470 

— 

— 

615 

63,266 

— 

— 

564 

106  ,420 

— 

— 

I 
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TABLE  17 

CRITICAL  SPEEDS  FOR  DAMAGED  WING 
AERODYNAMIC  DAMAGE  SITE  AT  90%  SEMI-SPAN 


No  Structural  Damage 


d i ve  rgen  ce 
speed 
mph 

D 

lbs 

f 1  utter 
speed 
mph 

D 

lbs 

1470 

0 

950 

0 

1029 

17,711 

1040 

18,090 

91  3 

27,887 

1850 

114,512 

776 

50,363 

— 

— 

699 

81  ,682 

— 

_ 

El -30%,  GJ-50% 

Decreased 

0  STA.5 ,6 

d i vergence 

D 

flutter 

r-s 

speed 

1  bs 

speed 

lbs 

mph 

mph 

950 

1  5,096 

940 

14,780 

725 

43,961 

— 

— 

667 

74,374 

— 

- . 
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TABLE  17  CONT'D 


El -40%,  GJ -60% 

Decreased 

@  STA.5,6 

di vergence 

D 

flutter 

D 

speed 

lbs 

speed 

1  bs 

mph 

mph 

1302 

0 

839 

0 

917 

14,065 

912 

13,912 

703 

41  ,333 

— 

— 

653 

71,285 

— 

— 

EI-50%,  G J - 70% 

Decreased  @  STA. 

5,6 

di vergence 

D 

f 1  utter 

D 

speed 

lbs 

speed 

lbs 

mph 

mph 

1226 

0 

789 

0 

872 

12,719 

864 

12,487 

674 

37,994 

— 

— 

634 

67,197 

— 

— 
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TABLE  17  CONT'D 

EI-70%,  GJ-70%  Decreased  @STA.5,6 


vergence 

speed 

mph 

D 

lbs 

flutter 

speed 

mph 

D 

lbs 

1226 

0 

789 

0 

857 

12 ,285 

864 

12,487 

659 

36,322 

_ 

_ 

621 


65,819 


TABLE  18 


FIRST  SIX  COUPLED  VIBRATION  MODES 
NO  STRUCTURAL  DAMAGE 


First  Mode 


1 . 

1  .  7861  E-2 

2. 

8. 1 4  2  4  E  —  2 

3. 

2 . 0927E-1 

4. 

4. 1 153E-1 

5. 

6.8321E-1 

6. 

1 .OOOOE+O 

Second 

Mode 

1  . 

-7.0912E-2 

2. 

-2.5918E-1 

3. 

-4.461 7E-1 

4. 

-3.9579E-1 

5. 

9.8291 E-2 

6. 

1  .OOOOE+O 

a>i  =  55.626  (rad/sec) 

8. 1218E-5 
2 . 0540E-4 
4. 3005E-4 
7.5596E-4 
1 . 0918E-3 
1 . 41 70E-3 

(i) 2  =184.114  (rad/sec) 

-2. 3049E-3 
-5 . 3055E-3 
-9 . 3808E-3 
-1 . 2708E-2 
-1 . 3433E-2 
-1  .  1  905F-2 


TABLE  18  CONT'D 


FIRST  SIX  COUPLED  VIBRATION  MODES 
NO  STRUCTURAL  DAMAGE 


Thi  rd 

Mode 

=  201.391  ( rad/sec) 

1  . 

1 . 1515E-1 

-1 . 3058E-2 

2. 

4 . 4077E-1 

-3.071 7E-2 

3. 

8. 3226E-1 

-5 . 5749E-2 

4. 

1 . OOOOE+O 

-8.6347E-2 

5. 

7 . 0826E- 1 

-1 . 1601E-1 

6. 

-1 . 5209E-2 

-1 .4203E-1 

Fourth 

Mode 

w4  =  355.122  (rad/sec) 

1 . 

5 . 37  31  E -  3 

2.1714E-2 

2. 

-3.0053E-2 

3 . 8287E- 3 

3. 

-1 . 3595E-1 

3. 3356E-2 

4. 

-1 . 9454E-1 

-5 . 9493E-3 

5. 

1 .2846E-1 

-5 .91 37E-2 

6. 

1 .OOOOE+O 

-1 .0816E-1 

Ik*  ... 
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TABLE  18  CONT'D 

FIRST  SIX  COUPLED  VIBRATION  MODES 
NO  STRUCTURAL  DAMAGE 


Fifth 

Mode 

u>5  =  429.720  (rad/sec) 

1 . 

-2.551 0  E - 1 

-3.2359E-3 

2. 

-6 . 5983E-1 

-1  . 2928E-3 

3. 

-3.9736E-1 

8 . 6078E-3 

4. 

6.5964E-1 

2.1318E-2 

5. 

1  .  OOOOE+O 

5 . 4350E-4 

6. 

-6. 3645E-1 

-6.5914E-2 

Sixth 

Mode 

u>6  =  542.953  (rad/sec) 

1  . 

6 . 0846E-2 

-2. 6348E-2 

2. 

1  . 9756E- 1 

-1 . 8475E-2 

3. 

5  .  7742E-2 

3.6184E-2 

4. 

-5 . 4663E-1 

5.014SE-2 

5.  ‘ 

-6. 3671 E-l 

-4.8803E-3 

6. 


1 .OOOO+O 


-7. 3730E-2 
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TABLE  19 


FIRST  SIX 

COUPLED  VIBRATIONAL  MODES 

El -30%, 

GJ -50%  DECREASE  @  STA.  5,6 

First 

Mode 

=  55.327  (rad/sec) 

1  . 

I . 7494E-2 

7.9823E-5 

2. 

7.9777E-2 

2.0194E-4 

3. 

2.0515E-1 

4 . 2304E-4 

4. 

4 . 0369E-1 

7. 4481E-4 

5. 

6. 73I5E-1 

1 .1935E-3 

6. 

1 .0000E+0 

1  .8490E-3 

Second 

Mode 

=  179.263  (rad/sec) 

1  . 

-7. 1562E-2 

-1  .6471  E-3 

2. 

-2 . 6485E-1 

-3. 7713E-3 

3. 

-4 . 6860E-1 

-6.6213E-3 

4. 

-4.5350E-] 

-8.5789E-3 

5. 

9. 1 557E-3 

-7.8630E-3 

6. 


1 . OOOOE+O 


-2. 1363E-3 
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TABLE  19  CONT'D 

FIRST  SIX  COUPLED  VIBRATIONAL  MODES 
El -30% ,  GJ -50%  DECREASE @  STA.  5,6 


Third  Mode 


w,  =  193.045  (rad/sec) 


1.  8 . 6587E-2 

2.  3.4417E-1 

3.  6 . 9724E-1 

4.  9.71 89E-1 

5.  1 . OOOOE+O 

6.  7 . 0264E- 1 


-1 .4718E-2 
-  3 . 4806E-2 
-6. 3928E-2 
-1 .001  IE-1 
-1 . 4886E-1 
-2.1 636E-1 


Fourth  Mode 


u>4=  320.403  (rad/sec) 


1.  -5 . 0885E-2 

2.  -1 . 9085E-1 

3.  -3.0635E-1 

4.  -2 . 1 647E-1 

5.  2.1238E-1 


1 . 7111E-2 
3. 3200E-2 
3 . 8754E-2 
2 . 1 283E-2 
-  3 . 0627E-2 


6. 


1 .0000E+0 


-1 . 3863E-1 


7o 


TABLE  19  CONT'D 

FIRST  SIX  COUPLED  VIBRATIONAL  MODES 
El- 30% ,  GJ -50%  DECREASE  @STA.5,6 


Fifth 

Mode 

=  388.648  (rad/sec) 

1  . 

-2.2937E-1 

-1 .0620E-2 

2. 

-6.1366E-1 

-1 . 4671 E-2 

3. 

-4 . 4284E- 1 

-2.21 36E-3 

4. 

5.1701E-1 

2.91 44E-2 

5. 

1  .  OOOOE+O 

4 . 3232E-2 

6. 

-7 . 7407E-1 

-5 . 5731 E-2 

Sixth 

Mode 

ug  =  503.248  (rad/sec) 

1  . 

1  .3173E-1 

-2.0362E- 2 

2. 

3. 5538E-1 

-2.1157E-2 

3. 

1  . 4274E- 1 

1 . 5050E-2 

4. 

-6 . 4707E-1 

4 . 6765E-2 

5. 

-8.9515E-1 

3.6104E-2 

6. 

1 .OOOOE+O 

-6.71 33E-2 
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TABLE  20 

FIRST  SIX  COUPLED  VIBRATIONAL  MODES 
EI-40%,  G J - 60%  DECREASE  @  STA.  5,6 


First 

Mode 

=  55.188  (rad/sec) 

1 . 

1 .7319E-2 

7.9137E-5 

2. 

7 . 8994E-2 

2 . 0024E-4 

3. 

2.0319E-1 

4. 1957E-4 

4. 

3 . 9996E- 1 

7. 3922E-4 

5. 

6.6816E-I 

1 . 2272E-3 

6. 

1 .OOOOE+O 

2 . 0493E-3 

Second 

Mode 

03 2  =  1  76 . 6  79  (  rad/sec) 

I  . 

-7 . 2892E-2 

-1 . 3041 E-3 

2. 

-2.71 57E-1 

-2 . 9632E-3 

3. 

-4.8737E-1 

-5.  1405E-3 

4. 

-4.9189E-1 

-6 . 2855E-3 

5. 

-4.4556E-2 

-4 . 2758E-3 

6. 

1 .OOOOE+O 

5.2893E-3 

TABLE  20  CONT'D 


FIRST  SIX  COUPLED  VIBRATIONAL  MODES 
EI-40%,  GJ -60%  DECREASE  0  STA.  5,6 


Third 

Mode 

oig  =  1  89 . 298  (  rad/sec) 

1  . 

4 . 9298E-2 

-1 . 2555E-2 

2. 

2.0447E-1 

-2.9741E-2 

3. 

4.4459E-1 

-5 . 4864E-2 

4. 

7.0332E-1 

-8 . 61 23E-2 

5. 

9 . 0890E- 1 

-1 . 3226E-1 

6. 

1 .OOOOE+O 

-2 . 0698E-1 

Fourth 

Mode 

=  305 . 886 ( rad/sec) 

1  . 

-6.4990E-2 

1 . 5677E-2 

2. 

-2 . 3483E-1 

3.  1  442E-2 

3. 

-3 . 6645E-1 

4 . 0023E-2 

4. 

-2 . 5528E-1 

3.001  IE-2 

5. 

2 . 0585E-1 

-1 . 5284E-2 

6. 


1 . OOOOE+O 


-2 . 5528E-1 


31 


TABLE  20  CONT ' D 

FIRST  SIX  COUPLED  VIBRATIONAL  MODES 
EI-40%,  G J - 60%  DECREASE  @STA.  5,6 


Fifth 

Mode 

tog  =  375 . 809  ( rad/sec) 

1  . 

— 2.1571E— 1 

-1 . 2203E-2 

2. 

-5 . 8362E-1 

-1 . 8276E-2 

3. 

-4.4140E-1 

-7.2774E-3 

4. 

4 . 7227E-1 

2 . 6343E-2 

5. 

1 . OOOOE+O 

5 . 5392E-2 

6. 

-8 . 4272E-1 

-4.4001E-2 

Sixth 

Mode 

oig  =  488 . 41  3(  rad/sec  ) 

1  . 

1  . 51 59E-1 

-1 .8832E-2 

2. 

4 . 0575E-1 

-2 . 1 748E-2 

3. 

1  .8497E-1 

8 . 8428E-3 

4. 

-6 . 6586E-1 

4.2804E-2 

5. 

-9 . 9924E-1 

5 . 0966E-2 

6. 

1  .OOOOE+O 

-6. 5087E-2 

TABLE  21 


be 


FIRST  SIX  COUPLED  VIBRATIONAL  MODES 
EI-50%,  GJ  -70%  DECREASE  (3  STA.  5,6 


First 

Mode 

a) i  -  55.006  (rad/sec) 

1  . 

1 . 7096E-2 

7 . 8320E-5 

2. 

7 . 7996E-2 

1  .9821  E-4 

3. 

2 . 0069E-1 

4 . 1  547E-4 

4. 

3.9520E-1 

7.3273E-4 

5. 

6.61 64E-1 

1 . 2742E-3 

6. 

1 .OOOOE+O 

2 . 3780E-3 

Second 

Mode 

oj2  =  1  7 2  .  7 5 3 (  rad/sec) 

1  . 

- 7 . 81  81 E -2 

-6 . 9945E-4 

2. 

-2 . 9435E-1 

-1  .  5207E-3 

3. 

5 . 4007E-1 

-2. 4379E-3 

4. 

-5. 7986E-1 

-1  .921  IE-3 

5. 

-1 . 487 1 E - 1 

3. 4304E-3 

6. 


1 . OOOOE+O 


2.4021E-2 


TABLE  21  CONT'D 


S3 


FIRST  SIX  COUPLED  VIBRATIONAL  MODES 
EI-50%,  GJ -  70%  DECREASE  OSTA.5,6 


Third 

Mode 

wg  =  183,658(rad/sec) 

1  . 

-5.3312E-4 

-6.9230E-3 

2. 

1 .0440E-2 

-1 . 6429E-2 

3. 

6 . 5776E-2 

-3.0463E-2 

4. 

2.1 757E-1 

-4. 7805E-2 

5. 

5 . 1 705E-1 

-7.6013E-2 

6. 

1 .OOOOE+O 

-1 .3159E-1 

Fourth 

Mode 

=  286.639(rad/sec) 

1  . 

-7.8687E-2 

1  .4291  E-2 

2. 

-2.81 12E-1 

2.9797E-2 

3. 

-4 . 4222E-1 

4.  1  693E-2 

4. 

-3. 3011E-1 

4 . 0233E-2 

5. 

1  . 6070E-1 

6 . 847 1  E -  3 

6. 


1 . OOOOE+O 


-1 . 4465E-1 
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TABLE  21  CONT'D 

FIRST  SIX  COUPLED  VIBRATIONAL  MODES 
El  -50%  ,  G  J  -  70%  DECREASE  l?  STA.  5,6 


F  i  f  th 

Mode 

cog  =  36 1 .  625  ( rad/ sec ) 

1 . 

-1 .9899E-1 

-1 . 3638E-2 

2. 

-5.4501E-1 

-2.1904E-2 

3. 

-4.3189E-1 

-1 . 3350E-2 

4. 

4 . 2563E-1 

2 . 0836E-2 

5. 

1  .  OOOOE  +  O 

6. 7331E-2 

6. 

-9.1 850E- 1 

-2.8823E-2 

Sixth 

Mode 

ug  =  469 . 745(rad/sec) 

1  . 

-1  . 5485E-1 

1 . 4987E-2 

2. 

-4. 1473E-1 

1 . 9435E-2 

3. 

-2 . 1  788E-1 

-1 . 529 1 E- 3 

4. 

6 . 0003E-1 

-3. 1527E-2 

5. 

1  .  OOOOE+O 

-5 . 9862E-2 

6, 

-8. 7734E-1 

5 . 3384E-2 
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TABLE  22 

The  Influence  of  Coupled  and  Uncoupled 
Vibrational  Modes  used  on  Critical  Speed 
for  Aerodynamic  Site  At  90%  Semi -Span 
with  0C  Sweep 
No  structural  damage 


di vergence 
speed 
mph 

f 1  utter 
speed 
mph 

A 

B 

A 

B 

1470 

1485 

950 

920 

1029 

•  1030 

1050 

965 

776 

777 

— - 

— 

699 

677 

— _ 

— 

El  30%,  GJ 

50%  Structural 

damage 

di vergence 
speed 
mph 

flutter 

speed 

mph 

A 

B 

A 

B 

950 

953 

940 

895 

725 

728 

— 

— 

667 

670 

— 

— 

NOTE:  A  Uncoupled  vibration  modes  used 

B  Coupled  vibration  modes  used. 


TABLE  22  CONT'D 


3.  El  40%,  GJ-60%  Structural  damage 


di vergence 
speed 
mph 

flutter 

speed 

mph 

A 

B 

A 

B 

1302 

1299 

839 

816 

917 

920 

912 

867 

703 

707 

— 

— 

653 

654 

— 

— 

E 1 50% ,  GJ  70% 

Structural  damage 

di  vergence 
speed 
mph 

flutter 

speed 

mph 

A 

B 

A 

B 

1226 

1242 

789 

770 

872 

874 

864 

820 

674 

679 

— 

— 

634 


636 


TABLE  22  CONT'D 


5.  El  70%,  GJ  70%  Structural  damage 


di vergence 
speed 
mph 

flutter 

speed 

mph 

A 

B 

A 

B 

1226 

1183 

780 

770 

857 

866 

864 

823 

657 

665 

— 

— 

621 

625 

_ 

- 

APPENDIX  A 


ESTIMATE  OF  DRAG  INCREMENTS  DUE  TO  DAMAGE 


A.l  Introduction 

Duri'-g  the  course  of  this  study  it  became  evident  that  the 
current  state  of  the  art  for  estimating  damage  induced  drag  levels 
leaves  much  to  be  desired.  For  the  estimation  of  the  influence  of 
aerodynamic  damage  on  the  performance,  the  drag  rise  should  be 
evaluated  on  the  basis  of  constant  lift  in  order  to  simulate  the  lift 
required  to  maintain  desired  flight  conditions.  This  basis  may  cause 
an  over-estimation  of  the  drag  in  the  present  study  of  aeroelastic 
failure.  Here,  possible  reduction  in  lift  should  also  be  considered 
in  the  analysis,  by  estimating  the  drag  rise  on  the  basis  of  fixed 
angle  of  attack.  Actual  damage  will  generally  require  aileron  deflec¬ 
tion  to  overcome  possible  asynwetries  in  the  rolling  moments  generated 
by  the  two  halves  of  the  wing.  This  asymmetry  depends  on  the  loss  in 
lift  due  to  the  aerodynamic  damage  and  its  location.  A  drag  increment 
based  on  fixed  angle  of  attack  takes  into  account  the  loss  in  lift 
and  is  believed  to  be  a  low  bound  since  it  does  not  include  the  lift 
associated  with  the  aileron  deflection.  The  two  estimates  for  ACg  due 
to  damage  are  shown  schematically  in  Figure  A.l.  It  is  believed  that 
the  actual  case  lies  between  estimates  based  on  fixed  lift  coefficients 
and  fixed  angle  of  attack.  A  better  estimate  cannot  be  made  at  the 
present  stage  of  the  study,  since  the  required  aerodynamic  data  is  lacking. 
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A.  2  Estimate  of  Drag  Increment  Due  to  Aerodynamic  Damage 

The  first  systematic  experimental  investigation  of  the  aero¬ 
dynamic  characteristics  of  damaged  wings  is  a  Cornell  Aeronautical 
Laboratory  report  (11)  dated  1952.  The  test  model  was  a  24"  chord  wing, 
with  NACA  65^012  airfoil,  that  spanned  the  102"  height  of  the  test 
section.  The  wing  was  tested  at  Mach  numbers  0.3,  0.7  and  0.85  with 
various  simulated  damage  configurations  that  are  shown  in  Figure  A. 2 
which  is  Table  I  of  Reference  11  . 

The  test  results  were  processed  in  a  later  report  by  the  same 
laboratory  (9)  by  fitting  a  parabolic  drag  model  to  the  data.  Using 
the  notation  of  this  reference: 


Co  =  CD  +  pCL’ 

IT 

changes  of  CQ  caused  by  damage  are  recorded  as  changes  of  the  "parasitic" 

2 

drag  coefficient,  C^,  and  the  slope  p  =  dCp/dC[_.  The  results  were  nor¬ 
malized  to  the  damage  hole  area,  instead  of  wing  area  by  introducing 

*  =  AD 
L0  "qA 

★  *  ? 

.  c0  tptL 

IT 

Figure  A. 3,  which  is  Table  II  of  Reference  9  ,  is  a  summary  of  drag 
increases  due  to  damage  holes  of  various  configurations.  This  data  is 
the  basis  for  the  estimations  made  in  Section  3,  4  and  5. 
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Before  we  proceed  and  discuss  the  results  of  this  investi¬ 
gation  it  is  important  to  make  notes  on  some  special  features  of  the 

test  conditions.  First,  the  basic  wing  is  two-dimensional,  thus  the 

★ 

slope  p*  in  the  expression  for  CQ  should  not  be  interpreted  as  induced 
drag  in  the  sense  of  a  finite  wing.  Secondly,  the  airfoil  had  a  thick¬ 
ness  ratio  of  12%  which  was  typical  of  aircraft  of  the  time  of  the  study 
This  value  is  much  higher  than  the  thickness  ratios  of  6%  or  less  used 
on  most  contemporary  fighters.  This  difference  is  of  importance  since 
transonic  phenomena  are  much  more  violent  for  the  thick  airfoil.  For 
example,  it  shows  a  dip  in  the  lift  curve  slope  at  a  Mach  number  about 
0.85,  an  early  drag  divergence  and  steeper  drag  rise  (12).  The  Mach  num¬ 
ber  0.85  tests  were  executed  only  over  a  very  narrow  range  of  angles  of 
attack,  namely  over  very  small  values  of  lift  coefficients.  Therefore, 
the  slope  term,  p  ,  was  not  evaluated  for  this  Mach  number.  The  discus¬ 
sion  below  will,  therefore,  be  limited  to  M  3  0.7.  The  investigation  in 
eluded  several  configurations  with  edges  of  the  holes  raised  to  form  a 
"scoop"  or  a  "spoiler"  lip.  The  height  of  the  lips  were  1/2"  and  1", 
or  about  2%  and  4%  of  the  local  chord.  Naturally,  such  protuberances 
generate  a  considerable  drag.  However,  according  to  information  re¬ 
ceived,  actual  petals  on  battle  damaged  wings  of  high  performance  fight¬ 
ers  do  not  exceed  1"  which  is  only  a  fraction  of  1%  of  full  size  typical 
chord.  Some  care  should,  therefore,  be  utilized  when  applying  damage 
cases  E  &  F  of  Figure  A. 2. 

The  damage  configurations  that  gave  the  highest  drag  rise 


are  1A15-2A15  and  1A10-2A10.  Recall  from  Figure  A. 2  that  hole  loca¬ 
tions  1  and  2  are  at  the  leading  edge  and  at  quarter  chord  upper  surface. 
It  is  believed  that  the  combination  of  a  leading  edge  hole  which  acts  as 
an  inlet,  a  hollow  wing  that  serves  as  a  settling  chamber  and  a  hole  lo¬ 
cated  in  the  region  of  maximum  suction  generate  a  fountain  which  acts  on 
the  external  flow  as  a  spoiler  that  causes  separation  and  the  associated 
drag  rise.  For  example,  for  damage  configuration  1A15-2A15  we  find 
A/S  =  0.00415  and 


C  =24 

TT 


p  =  10.0 


which  give  for  a  fixed  Cp  =  0.8  a  drag  rise  of  Cp  =  8.8  or  its  equivalent 
ACp  =  0.036.  The  alternative  estimate,  namely  the  one  bated  on  fixed 
angle  of  attack,  predicts  increase  in  drag  coefficient  is  0.017  which 
is  less  than  half  the  value  predicted  for  constant  CL> 

A  through  hole  is  represented  by  damage  configuration  2A15- 
4A15  whose  two  holes  are  centered  at  quarter  chord.  For  this  configur- 


a  ti  on 


CD*  -  0 

7T 


★ 

p  =  6 


For  =  0.8  the  predicted  increase  in  drag  coefficient  is 

★ 

Cp  =  3.84  or  ACp  =  0.016.  For  the  case  of  constant  angle  of  attack  we 
find  for  the  damaged  wing  at  a  =  4.6°,  =  0.64  and  CQ  =  0.028,  i.e. 

an  increase  of  ACp  =  0.006  v/h i ch  again  is  less  than  half  the  value 


To  conclude  the  analysis  of  Reference  11  and  9,  we  summarize 


findings  for  CL  =  0.8  (undamaged)  and  M  =  0.7. 


CONFIGURATION  REPRESENTING  flxed  CL  flxeda  fixed  CL  fixed  a 


1A15-2A15 

Leading  edge  inlet 
forms  a  fountain 

8.8 

4.1 

.036 

.017 

2  A 1 5  -  4  A 1 5 

Through  hole'at 
quarter  chord 

3.8 

1.5 

.016 

.006 

However,  we  recall  that  these  values  are  not  necessarily  valid  in 
transonic  speeds,  that  they  do  not  include  the  expected  large  change 
in  induced  drag  which  results  from  the  modification  in  lift  distribu¬ 
tion  and  that  they  do  not  include  the  additional  effects  associated 
with  aileron  deflection  that  should  follow  any  non-symmetrical  damage. 
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Table  I 

* 

Key  to  Configuration  Notation 


Hole  Location 

Hole  Number  Location 

1  Leading  edge 

2  .25  chord,  upper  surface 

3  .70  chord,  upper  surface 

4  .25  chord,  lower  surface 

5  .70  chord,  lower  surface 


Lip  Configuration 

Symbol  Type  of  Lip  Depth  of  Flange 


A 

flush 

none 

B 

protruding 

1/2 

in. 

C 

receding 

1/2 

in. 

D 

scoop 

1/2 

in. 

E 

scoop 

1 

in. 

F 

spoiler 

1/2 

in. 

G 

spoiler 

1 

in. 

Description  of  Hole  Size 


Example  of  Various  Lips 


Lip  "A"  for  a 
loading  edge  hole 


stream 


Lip  "A"  for  other  than 
leading  edge  holes 


Numeral 


Hole  Diameter* 


10 

15 


.10  chord 
.15  chord 


*For  leading  edge  holes,  the  diameter 
of  an  equal  area  circular  hole. 


Lips  "F"  &  "G" 


stream 


%  on  Lip  "F" 
1  on  Lip  ”G” 


denotes  wing  interior 


DAMAGE  CONFIGURATIONS  TESTED  IN  REFERENCE  11 
Figure  A. 2 


SUMMARY  OF  DRAG  INCREASES  ATTRIBUTABLE  TO  DAMAGE  HOLES  IN  WINGS 


APPENDIX  B 


PROBLEM  FORMULATION 


The  classical  thin  beam  buckling  equations  employed  in  the 
present  analyses  were  first  developed  by  Michel  1  and  are  currently 
derived  in  many  classical  texts  on  structural  analysis  (13).  The  math 
model  utilized  for  the  present  study  is  illustrated  in  Figure  4.  Here 
the  reference  axes  are  taken  along  the  elastic  axis  of  the  wing  which 
may  be  at  an  angle  of  sweep  A  and  a  root  chord  angle  of  attack  aQ 
relative  to  the  on-coming  flow.  In  addition,  the  damage  is  assumed  to 
be  localized  giving  rise  to  modifications  in  the  aerodynamic  loads  that 
can  be  represented  as  localized  changes  in  drag,  lift,  and  moment.  These 
concentrated  changes,  illustrated  in  Figure  2,  are  referenced  to 
chordwise  coordinate  n  forward  or  aft  of  the  elastic  axis  and  spanwise 
coordinate  £.  Writing  the  equilibrium  equations  in  the  deformed  state 
of  the  wing  yields: 

EI)1  ^7  =  M r  *  El),  ^7 
1  dr  ^  1  dxc 

El)  d4  =  -  Mn  «  EI)?  (IB) 

*  dr  n  d*x 


GJ) 


eff 


d0 


=  M, 


*  GJ)eff 


For  fighter  wing  configurations  the  chordwise  stiffness  is  much  greater 
than  the  stiffness  normal  to  the  chord  direction,  i.e. 
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w*»- 
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EI)l 

eT5^ 


»  1 


and  the  last  two  equations  essentially  uncouple,  to  the  first  order,  from 
the  first  one.  This  is  due  to  the  smallness  of  the  chordwise  deflection 
n  and  its  lower  order  derivatives  when  compared  to  the  wing  twist  0 
and  deflection  £  normal  to  the  wing  chord  plane.  That  is,  for  the  small 
angles  and  deflections  occurring  at  the  onset  of  a  flutter  or  divergence 
instabil ity 

n  *  v  c  »  to  0  -  b 

with  3  projection  of  0  on  y-z  plane 
v  projection  of  n  on  y-z  plane 
to  projection  of  <;  on  y-z  plane 


] 

and  deflections  along  the  x  axis  are  neglected.  Further  simplifications 
can  be  made  upon  assuming  curvatures  in  n£  and  planes  are  very 
nearly  equal  to  the  curvatures  in  the  xy  and  zx  planes  respectively 
and  also  that 


-d!  = 
d£  dx 
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Displacements  are  also  considered  small  (flutter  or  divergence  just 
initiated)  cosines  of  all  angles  are  equal  to  unity  and  sines  of  all 
angles  are  equal  to  the  angle  and  its  tangent.  Bending  moments  and 
torques  in  the  cross  section  of  the  deformed  configuration  can  be 
referenced  to  the  original  coordinate  state  by  the  following  transformation 


Now  the  moments  M  ,  M  ,  and  M  created  by  the  damage  induce  aerodynamic 
x  y  z 

forces  are  from  Figure  4. 

Mx  =  -  AD  cos  A  [u>(0  +  nB(0  -  w(x)] 

-  AL  [v(0  +  n  -  v(x)  -  7]  -—]  +  AM 

My  =  AD  sin  A  [(40  +  oB(0  -  w(x)]  (3B) 

+  al  [U  -  x)  -  n~  (03 

Mz  =  AD  cos  A  [U  -  x)  -  (C) n  1 


-  AD  sin  A  [v(0  +  n  -  v(x)] 


Taking  the  appropriate  derivatives  of  the  second  and  third  equations  of 
IB  provides,  along  with  2B: 


El 


d2m 

s? 


+ 


U  H  it. 

7 r +  7?  [3MZ] 

dx  dx 


=  L 


(4B) 


d_ 

dx 


W  +  37 


z  ax 


=  M 


where  L  and  M  are  the  distributed  lift  and  moment  on  the  lifting  surface. 

Coupling  (3B)  and  (4B)  gives  rise  to  damage  induced  lift  and  moment 
terms  dropping  out  as  higher  order  effects,  while  only  the  damaged  induced 
drag  terms  are  remaining.  This  result  is  also  based  upon  the  obser¬ 
vations  from  wind  tunnel  tests  that  AL  and  AD  are  roughly  of  the 
same  order.  At  large  angles  of  attack  &0  the  change  in  lift  term  AL, 
which  is  usually  negative,  may  contribute  a  small  amount  to  the  chord- 
wise  drag  force  tending  to  diverge  the  wing.  For  moderate  to  small  angles 
of  attack,  however,  this  effect  is  small  and  the  equations  take  the 
form: 


Q 

dx 


+  AD  cos 


at 


at 


2 

A  (e-x)  -5-7  -  AD  sin  A  n  ^  =  M 
ax  ax 


dx 


El 


d2o 

d? 


at  a  at* 


-  AD  sin  A  4-7  -  2AD  cos  A 

dx  2  2 

+  AD  cos  A  (£-x)  --4  -  AD  An  sin  -d— f  =  L 
dx*  dx* 


The  above  constitutes  the  basic  equations  employed  for  the  flutter 


and  divergence  analyses  of  the  present  study. 


APPENDIX  C 


THREE-DIMENSIONAL  FINITE  ELEMENT  MODELING  QF  STRUCTURAL 
DISCONTINUITIES  IN  CONSTANT  CHORD  MINGS 
OR  ROTOR  ELEMENTS 


To  illustrate  the  detail  in  finite  element  structural  modeling  that 
can  be  accomplished  for  a  well  defined  structure,  reference  is  made  here 
to  an  earlier  investigation  of  structural  discontinuities  on  constant  chord 
rotors  or  wings  with  a  single  cell  torque  box.  The  study  is  also  an 
alternate  check  on  the  influence  that  localized  structural  damage  has  on 
the  mode  shape  frequency  characteristics  of  undamaged  wings  or  rotors. 

A  three-dimensional  finite  element  model  of  a  typical  composite 
rotor  blade  has  been  developed  in  Reference  20  to  investigate  the  influence 
of  structural  repair  discontinuities  reflected  in  the  blade's  stiffness 
and  model  characteristics.  Modern  computational  procedures  including 
multi-level  substructuring  (Reference  21)  are  employed  during  the  solution 
of  symmetric  finite  element  equations.  The  computed  stiffness  and  modal 
characteristics  of  the  blade  are  correlated  to  those  obtained  from  labora¬ 
tory  tests.  Also  included  in  this  report  are  the  results  of  the  pre¬ 
liminary  investigations  of  the  influence  of  idealized  repair  discontinuities 
as  reflected  by  the  finite  element  model. 

As  can  be  seen  in  Figure  Cl-a,  the  rotor  blade  is  fairly  uniform 
in  cross-section  along  its  entire  length  of  22  ft.  with  the  exception  of 
the  additional  grip  plates  and  doublers  overlaid  on  the  4  ft.  portion  at 
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the  supported  end.  During  the  laboratory  tests,  the  blade  was  mounted  on 
a  relatively  rigid  support  fixture. 

In  the  finite  element  model,  the  blade  was  assumed  to  have  constant 
cross-sectional  properties  within  two  spanwise  regions  as  shown  in  Fig. 

Cl-b.  The  difference  between  these  regions  is  due  to  the  fact  that  thicker 
membrane  elements  were  required  near  the  support  so  that  the  effect  of  the 
grip  plates  and  doublers  (Fig.  Cl-a)  could  be  reflected  in  the  idealization. 
As  shown  in  Fig.  Cl-c,  three-dimensional  brick  elements  were  used  to  model 
the  spar  and  the  honeycomb  core  while  membrane  elements  were  used  for  the 
skin.  The  typical  section  in  each  spanwise  region  has  21  elements  (9  bricks 
and  12  membrane  elements)  and  44  nodal  points  (Fig.  Cl-c).  The  total  finite 
element  idealization  (Fig.  22)  consists  of  88  of  these  typical  sections  each 
being  assumed  to  have  identical  stiffness  and  mass  properties  within  each 
spanwise  region.  This  results  in  1848  elements  and  1958  nodal  points  with 
three  degrees-of- freedom  at  each  node. 

The  reduction  process,  which  is  a  simultaneous  static  and  Guyan  reduc¬ 
tion  of  stiffness  and  mass  matrices,  respectively,  takes  advantage  of  the 
repeated  use  of  identical  substructures  (Ref.  21).  First  the  21  elements 
of  Fig.  Cl-c  were  processed  in  each  spanwise  region  to  produce  two  of  the 
88  substructures  of  Fig.  3.  Subsequently  i'wo  of  these  substructures  were 
combined  in  each  region  to  generate  two  of  the  44  substructures  of  Fig.  C4. 
This  proces's  was  applied  two  more  times  to  produce  the  substructures  of 
Fig.  C5  and  C6.  The  11  elements  of  Fig.  C6  were  then  combined  to  generate 
the  final  substructure  having  30  nodal  points  as  shown  in  Fig.  C7. 
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The  stiffness  and  mass  matrices  of  the  final  substructures  (Fig. 

C7)  were  then  used  to  determine  the  flexibility  and  modal  characteristics 
of  the  rotor  blade.  Fig.  C8  depicts  the  correlation  of  the  measured 
(Ref.  22)  and  computed  flexibility  coefficients  corresponding  to  the  verti¬ 
cal  deflection  of  the  blade  under  the  applied  unit  loads  as  shown  in  the 
figure.  The  discrepancy  was  less  than  1.5  percent.  Table  C-l  gives  the 
correlation  of  the  measured  (Ref.  23)  and  computed  natural  frequencies 
of  the  blade.  Also  included  in  this  table  are  the  computed  frequencies 
from  a  beam  model  (Ref.  23)  that  requires  predetermined  cross-sectional, 
stiffness,  and  mass  characteristics  of  the  blade.  Although  the  finite 
element  model  gave  larger  discrepancies  than  the  beam  model  for  the  out- 
of-plane  bending  modes  in  this  particular  case,  it  could  be  found  attrac¬ 
tive  for  the  geometric  generality  it  provides  in  the  idealization  of  local 
structural  discontinuities  such  as  holes,  repairs,  stiffeners,  etc. 

Currently,  the  influences  of  the  three  types  of  damage  shown  in 
Fig.  C-9  on  the  stiffness  and  modal  characteristics  of  the  blade  have  been 
investigated  as  a  preliminary  study  of  more  general  repair  discontinuities. 
The  first  two  damage  types,  i.e.  HOLE  and  NOTCH,  result  from  complete 
delamination  of  the  core  and  skin  on  both  faces  of  the  blade  adjacent  to 
the  leading  edge  torque  box  as  shown  in  Fig.  C9~a  and  C9~b.  In  this  study, 
the  spanwise  location  for  these  damages  has  been  selected  to  be  at  STA  75 
or  213  inches  from  the  tip.  The  third  damage,  i.e.  CRACK,  is  formed  by 
the  delamination  of  the  skin  on  both  faces  and  located  between  STA  72  and 
STA  264  as  shown  in  Fig.  C9-c.  This  damage  idealizes  the  separation  of  the 
leading  edge  torque  box  from  the  trailing  edge  torque  box  as  a  result  of 
longitudinal  cracking  of  the  skin  on  both  faces  of  the  blade. 
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The  effect  of  the  above  types  of  damage  on  the  flexibility  of  the 
blade  is  shown  in  Fig.  CIO.  The  load-deflection  figures  shown  in  this 
figure,  which  could  be  considered  a  computer  simulation  of  the  droop  test, 
correspond  to  computed  tip  deflections  of  the  blade  under  upward  applied 
forces  at  STA  288.  The  deflection  values  and  percent  changes  in  the 
deflection  due  to  the  damage  are  given  in  Table  C-2  for  all  stations 
considered  in  the  standard  droop  test. 

The  computed  natural  frequency  response  of  the  damaged  blades 
are  given  in  Table  C-3  for  the  first  five  vibration  modes.  Also  included 
in  this  table  is  the  percent  change  in  the  vibration  frequencies  due  to 
the  damage. 


retantion 


Finite  Element  Idealization  "f  the  Rotor  Blade. 
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Region  1.  Region  2. 

|* —  4  ft — *j-« - - - — 


18  ft 


Fig.  C-2.  Plan  View  Showing  the  Two  Typical  Cross- Sect  ions . 
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Fig.  C-3.  Reduction  to  88  Substructures. 
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Fig.  C-4.  Reduction  to  44  Substructures. 
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Fig.C-5.  Reduction  to  22  Substructures. 
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Fig.  C-6.  Reduction  to  11  Subs  tructures . 
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Fig.  C-7.  Reduction  to  Final  Substructure  with  30  Xodes 
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Fig.  C-8.  Comparison 
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BLADE  MODE 


First  Out-of-Plane 


First  In-Plane 
Second  Out-of-Plane 
Third  Out-of-Plane 
First  Torsion 


Experimental  Computed  Freq  &  %  Discrepancy 
Frequency 

DF1758  (Ref. 3)  Finite  Element 
•  HZ  HZ  7.  HZ  7. 


1.37 


DF1758  (Ref. 3) 

!  m 

1.44  j 

6.32 

5.3 

8.05 

0.6 

20.88 

4.9 

26.45  1.7  26.02 


Table  C-l.  -  Correlation  of  Calculated  and  lie a sued  Natural 
Frequencies . 
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Computed  Deflection  and  7.  Change  in  Deflection 
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(a)  HOLE  :  Delamination  of  core  and  skin  on  both  faces. 
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(b)  NOTCH  :  Delamination  of  core  and  skin  on  both  faces. 
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(c)  CRACK  :  Delamination  of  skin  on  both  faces. 


Fig.  C-9.-  Idealized  Damage  Types. 


